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SUMMARY 


A  series  of  flight  trials  was  performed  on  the  F-111C  aircraft  at  the  RAAF’s  Aircraft 
Research  and  Development  Unit  in  February  and  October  1987.  Data  obtained  from  the 
tests  were  analysed  at  the  Aeronautical  Research  Laboratory  to  determine  the  aircraft 
aerodynamic  and  control  derivatives.  This  report  describes  the  methods  and  computer 
programs  which  are  used  to  process  and  analyse  the  flight  test  data.  Data  handling 
procedures,  pre-analysis  flight  data  processing  and  the  methods  used  to  make  corrections 
to  air  sensor  measurements  are  described.  Although  the  test  programme  was  conducted 
on  a  F-111C  aircraft,  with  minor  alterations  the  computer  programs  and  procedures  can 
be  used  for  other  aircraft  test  programmes. 
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Introduction 


Models  of  aircraft  flight  dynamics  are  used  for  a  range  of  applications,  for  example, 
the  analysis  of  aircraft  behaviour,  for  aircraft  design,  and  for  driving  flight  simula¬ 
tors.  Part  of  this  modelling  procedure  involves  the  estimation  of  the  aerodynamic 
stability  and  control  derivatives  from  various  sources  such  as  empirical,  wir.d  tun¬ 
nel  and  flight  tests.  This  report  describes  the  methods  and  computer  programs 
which  are  used  to  process  and  analyse  data  obtained  from  a  flight  test  programme 
conducted  on  a  General  Dynamics  F-111C  aircraft.  Data  handling  procedures,  pre¬ 
analysis  flight  data  processing  and  the  methods  used  to  make  corrections  to  air 
sensor  measurements  are  described.  Data  extraction  i3  carried  out  on  a  VAX  750 
computer  using  Fortran  programs  developed  at  the  Aircraft  Research  and  Develop¬ 
ment  Unit  (ARDU).  Data  processing  and  analysis  is  carried  out  on  the  ARL  EIXSI 
6400  computer  using  Fortran  programs  either  acquired  or  specifically  developed  at 
ARL.  Data  presentation  is  carried  out  on  an  IBM  PS2  personal  computer  using  a 
computer  program  written  in  Pascal. 

Flight  testing  of  the  aircraft  was  carried  out  at  ARDU  in  1987.  Data  from  the 
programme  has  been  used  to  validate  a  comprehensive  flight  dynamic  model  of  the 
F-lllC  which  was  developed  at  ARL  and  this  model  will  be  used  to  upgrade  the 
F-111C  flight  training  simulator  at  RAAF  Base  Amberley. 


2  Test  Aircraft  and  Instrumentation 

The  F-lllC  test  aircraft  A8-132,  illustrated  in  Figure  1,  was  extensively  modified 
under  ARDU  Test  Schedule  1650  with  flight  test  quality  instrumentation  and  data 
recording  equipment.  This  equipment  known  as  the  Airborne  Flight- Test  Recording 
and  Analysis  System  (AFTRAS)  provides  on-board  digital  magnetic-tape  recording 
and  telemetry  information  for  real-time  flight  test  monitoring.  The  instrumentation 
was  developed  for  use  for  store-carriage  and  release  tests  and  for  the  flight  dynamic 
measurements  and  is  capable  of  measuring  200  measurands  at  a  sampling  rate  of 
60  per  second.  Special  equipment  was  developed  for  recording,  and  for  manual  ad¬ 
justment  of  the  pitch  and  roll  adaptive  gain  values.  Because  of  insufficient  time, 
the  instrumentation  required  to  monitor  other  parts  of  the  adaptive  control  system 
was  not  installed.  Similarly,  the  instrumentation  for  monitoring  engine  parame¬ 
ters  required  for  detailed  performance  measurements  was  not  fitted.  These  parts 
of  the  instrumentation  were  not  required  for  the  estimation  of  the  aerodynamic 
coefficients,  but  were  required  for  the  related  investigation  of  the  adaptive  control 
system  behaviour  and  aircraft  performance  characteristics. 

The  Nose  Boom  Transducing  Unit  which  was  used  for  the  second  phase  of  the 
trials,  was  constructed  by  the  Advanced  Engineering  Laboratory  (AEL)  Salisbury 
and  was  designed  to  provide  high  quai’ty  measurements  of  pitot  pressure,  angle- 
of-attack,  angle-of-sideslip  and  linear  accelerations  parallel  and  normal  to  the  local 
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airflow  direction.  The  unit  was  modelled  on  the  CONRAC  Nose  Boom  Instrumen¬ 
tation  Unit  (NBTU)  Model  25126F,  developed  by  the  USAF  for  flight  dynamic 
performance  measurements.  The  NBTU  is  shown  in  Figure  2  and  a  detailed  de¬ 
scription  of  the  assembly  is  given  in  Reference  [lj. 

A  list  of  the  instrumentation  channels  (measurands),  including  their  ranges  and 
accuracies,  which  are  used  for  the  flight  dynamic  analysis  is  given  in  Table  1.  These 
channels  were  recorded  at  a  rate  of  60  samples  per  second  using  the  AFTRAS 
system. 

The  measurement  of  the  pressure  error  corrections  for  the  NBTU  was  carried 
out  by  ARDU  as  part  of  the  Test  Schedule  1691.  Details  of  these  measurements  are 
presented  in  Appendix  D. 


3  Flight  Test  Programme 

An  outline  of  the  flight  test  programme  is  prer  ented  below.  Full  details  of  the 
programme  are  contained  in  the  flight  data  analysis  reports  (References  [2]  to  (7j). 

3.1  Test  Points 

The  flight  tes>.  programme  was  carried  out  in  two  phases.  The  matrix  of  test  points 
covered  in  Phase  i  .ad  Phase  2  is  given  in  Table  2.  The  first  phase  revered  75  test 
conditions  representing  combinations  of  wing  sweep,  altitude  and  Marh  number.  At 
each  test  point  two  longitudinal  and  two  lateral  manoeuvres  were  performed.  The 
flow  angles  in  this  phase  were  obtained  from  tranducers  in  the  aircraft’s  standard 
Central  .Air  Data  System  (CADS).  For  the  second  phase  of  the  programme  a  Nose 
Boom  Transducing  Unit  (NBTU)  was  available  for  the  measurement  of  the  angles 
of  attack  and  sideslip.  A  number  of  Phase  1  test  point  manoeuvres  were  repeated 
in  Phase  2  to  compare  the  accuracy  of  the  results  using  the  aircraft  system  (CADS) 
and  the  NBTU  system  for  measurements  of  angle  of  attack  and  angle  of  sideslip. 
The  remainmg  tests  were  made  at  Mach  numbers  between  those  tested  in  Phase  1 
to  provide  a  more  comprehensive  data  coverage.  The  total  test  program  required 
24.6  test  hours  of  flying. 

3.2  Teat  Manoeuvres 

The  following  manoeuvres  were  performed  at  each  test  point. 

1.  accurate  trim 

2.  pitch  input  where  the  stick  is  pulled  bad  (1  to  «  2  ’g’)  then  pushed 
forward  to  n<  utral  and  the  aircraft  allowed  to  damp  in  pitch. 

3.  trim 
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4.  pitch  input  where  the  stick  is  pushed  forward  (1  to  «  0  ’g’)  then  pulled 
back  to  neutral  and  the  aircraft  allowed  to  damp  in  pitch. 

5.  accurate  trim 

6.  rudder  step  input  to  left  followed  by  aileron  doublet  to  achieve  «  x30° 
bank  angle.  Rudder  and  aileron  released  together. 

7.  trim 

8.  manoeuvre  6.  repeated  but  rudder  input  to  right  and  opposite  roll  ap¬ 
plied 

9.  trim 

These  specified  manoeuvres  are  designed  to  give  an  aircraft  response  which  is  opti¬ 
mum  for  the  determination  of  stability  and  control  derivatives  using  the  techniques 
described  in  Section  4.  Advice  provided  by  NASA  Dryden  test  personnel  from  ex¬ 
perience  with  the  F-111A  TACT  aircraft,  indicated  that  large  rapid  control  inputs 
were  necessary  to  provide  large  amplitude  excitation  of  the  natural  modes  of  the 
aircraft  with  the  automatic  flight  control  system  engaged. 

Manoeuvres  were  also  flown  with  the  aircraft  in  the  landing  and  take  oil  config¬ 
urations.  These  cases  are  summarised  in  Table  3. 

Supplementary  manoeuvres  were  performed  to  enhance  the  prediction  capability 
of  the  validated  ARL  flight  dynamic  model  and  for  use  as  test  manoeuvres  for  the 
F-1I1C  flight  training  simulator.  These  included: 

1.  longitudinal  roller-coaster  manoeuvres 

2.  lateral  oscillatory  manoeuvres 

3.  dutch  rolls 

4.  steady  heading  sideslips 

5.  longitudinal  trims 

The  manoeuvres  and  the  flight  conditions  at  which  they  were  performed  are  sum¬ 
marised  in  Table  4. 

3.3  Flight  Control  System  Status 

The  tests  were  conducted  with  the  flight  control  system  in  the  normal  mode  and  with 
the  system  gains  determined  by  the  normal  adaptive  mode  gain  changer.  However, 
when  rapid  control  inputs  were  applied  at  some  flight  conditions,  motion  due  to  the 
adaptive  mode  was  superimposed  on  the  natural  motion  of  the  aircraft.  The  problem 
of  adaptive  mode  ringing  was  overcome  by  pumping  the  stick  and  driving  the  control 
system  gains  down  to  a  level  where  the  aircraft  natural  response  dominated  the 
motion. 
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Flight  Data  Processing  and  Analysis 


Flight  data  processing  and  analysis  was  carried  out  in  the  order  shown  in  Figure 
3.  The  procedures  used,  and  the  software  developed  and  acquired  for  this  purpose 
are  summarised  in  this  section.  Details  of  the  procedures  are  given  in  the  Appen¬ 
dices.  Figures  4  and  5  list  the  names  of  the  input  and  output  files  used  by  various 
procedures  and  lists  the  analysis  programs. 

A  significant  aspect  of  the  programme  involved  the  managment  of  large  amounts 
of  data.  Approximatley  10.5  million  data  points  were  acquired  in  this  programme, 
being  the  aggregate  of  33  measurands  sampled  at  60  samples  per  second  for  test 
manoeuvres  totaling  46  seconds  at  each  of  the  162  test  conditions.  Approximately 
20  hours  of  CPU  time  and  15  Megabytes  of  disk  storage  on  a  high  speed  mainframe 
computer  was  required  to  analyse  the  cases  for  a  complete  Mach  number  range  at 
one  sweep  angle  and  one  altitude.  A  total  of  6  sweep  angles  was  tested  at  5  different 
altitudes  in  the  clean  configuration.  Additional  measurements  were  also  made  in 
the  take  off  and  landing  configurations. 

4.1  Data  Extraction 

The  AFTRAS  flight  data  system  provides  organised  procedures  for  accessing  se¬ 
lected  channels,  for  applying  calibrations  to  give  engineering  units,  and  for  for¬ 
matting  the  data  for  subsequent  analysis.  Calibrations  are  stored  as  polynomial 
coefficients  against  the  date  of  calibrations  and  are  appended  to  the  data  files  ob¬ 
tained  from  each  flight.  Details  of  these  procedures  are  given  in  Reference  [8].  The 
procedures  are  carried  out  using  program  EXTRACT  which  can  be  run  on  the 
ARDU  or  ARL  VAX  computers.  Appendix  A  describes  the  EXTRACT  procedures 
in  detail  and  through  the  use  of  an  example,  indicates  the  information  which  has 
to  be  provided  and  the  form  of  the  data  extracted. 

The  input  data  required  for  use  by  the  parameter  estimation  techniques  needs 
additional  processing  to  provide  the  required  data  format  and  also  to  apply  further 
measurement  error  corrections.  This  processing  stage  is  carried  out  within  program 
FDP  on  the  ARL  ELXSI  computer  and  is  described  in  detail  in  Appendix  C.  The 
corrections  are  required  : 

1.  to  compensate  the  airspeed,  Mach  number  and  altitude  for  pressure  er¬ 
rors  and  compressibility  effects  at  the  pitot-static  sensor  locations 

2.  to  apply  time  shifts  to  the  time  history  records  to  compensate  for  instru¬ 
ment  signal  conditioning  and  recording  lags 

3.  to  calculate  pitch  and  roll  control  deflections  from  the  measured  stabila- 
tor  deflections 

4.  to  calculate  weight,  c.g.  and  moments  of  inertia  information  from  the 
fuel  tank  contents  data 
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4.2  Aircraft  Mass  Characteristics 


Accurate  information  is  required  of  all-up  weight,  horizontal,  vertical  and  lateral 
centre  of  gravity  position  and  inertia  properties  for  the  analysis  of  dynamic  ma¬ 
noeuvres.  This  information  is  used  to  convert  the  aerodynamic  derivatives  into 
non-dimensional  coefficients,  to  provide  references  for  instrumentation  and  sensor 
position  and  to  allow  vaiid  comparisons  of  the  derivatives  obtained  from  different 
flights  where  the  mass  characteristics,  especially  the  C  of  G  position,  may  be  dif¬ 
ferent.  Aerodynamic  moment  data  obtained  from  wind  tunnel  tests  are  usually 
referenced  to  a  given  C  of  G  position  and  an  adjustment  for  the  actual  flight  test 
positions  must  be  included  if  a  valid  comparison  is  to  be  made. 

Procedures  for  determining  the  mass  characteristics  for  the  test  aircraft  A8-132 
were  as  follows: 

•  Aircraft  was  weighed  to  establish  baseline  data. 

•  Fuel  calibration  tests  were  conducted. 

•  Software  was  developed  using  manufacturers  data  to  calculate  mass  charac¬ 
teristics  and  adjustments  made  to  reflect  weight  and  fuel  calibration  results. 

•  At  each  test  point  the  indicated  fuel  tank  contents  were  recorded. 

4.2.1  Aircraft  Weighings 

The  aircraft  was  weighed  in  accordance  with  standard  RAAF  weighing  procedures. 
Tables  5  to  10  show  the  results  for  the  various  configurations  which  form  the  baseline 
aircraft  information.  The  change  in  weight  and  C  of  G  for  the  aircraft  with  the 
NBTU  fitted  (Phase  2  flights)  was  derived  from  the  information  contained  in  Table 
10. 


4.2.2  Fuel  Calibration 

A  fuel  calibration  was  obtained  by  emptying  the  fuel  tanks  of  all  useable  fuel  and 
adding  a  known  amount  to  the  forward  and  aft  tanks  then  recording  the  indicated 
and  actual  tank  contents.  During  the  flight  test  the  aircraft  loading  and  centre-of- 
gravity  varied  in  accordance  with  the  aircraft  auto-fuel  schedule,  which  is  described 
in  Appendix  B.  After  take-off  the  aircraft  would  use  all  fuel  in  the  wing  and  primary 
weapons  bay  fuel  tan  Vs  before  any  flight  test  manoeuvres  commenced.  Prior  to  each 
test  point  the  indicated  foward  and  aft  fuel  tank  contents  were  recorded  on  the  pilots 
test  card  and,  by  applying  the  fuel  calibration,  an  accurate  calculation  of  C  of  G 
and  AUW  was  made.  Figure  6  shows  a  typical  test  card  with  relevant  information 
highlighted.  During  the  test  manoeuvres,  the  aircraft  weight  was  typically  near 
70,000  lb. 


4.3  Instrumentation  Time  Lags 

A  procedure  for  identifying  the  reiative  time  lags  between  instrumentation  channels 
was  developed  at  ARL  and  is  documented  in  Reference  [9j.  This  procedure  also 
uses  the  maximum  likelihood  technique  and  was  applied  to  a  number  of  selected 
time  histories  to  determine  the  lag  parameters.  The  resulting  parameters  are  given 
in  Table  11  where  a  positive  integer  indicates  n  time  samples  (n/60  seconds)  lagged 
with  respect  to  the  control  deflection  signals.  The  table  shows  that  the  NBTU 
signal  for  angle-of-attack  leads  the  control  deflection  signals  by  two  sample  inter¬ 
vals,  showing  that  this  instrumentation  has  smaller  signal  delays  than  the  Phase  1 
instrumentation. 

4.4  Calibration  of  a  and  / 3  Flow  Vanes  Using  Flight  Path 
Reconstruction 

In  addition  to  determining  the  pressure  error  corrections  to  airspeed  and  altitude  it 
is  also  necessary  to  determine  the  position  errors  (or  calibration  constants)  for  the 
angle-of-attack  and  angle-of-sideslip  tranducers.  In  particular  the  aircraft  CADS 
transducers  are  mounted  close  to  the  forward  fuselage  where  local  flow  angles  can 
differ  substantially  from  the  free-stream  value.  A  Flight  Path  Reconstruction  (FPR) 
method  for  determining  these  calibration  constants  is  described  in  Reference  [10] . 

The  flight  path  reconstruction  technique  uses  an  extended  Kalman  filter  to  de¬ 
termine  the  calibration  parameters  relating  the  computed  and  measured  output 
variables. 

This  method  uses  a  combined  parameter  and  state  estimation  technique  and 
is  based  on  using  the  kinematic  equations  of  motion  of  a  rigid  body.  Using  longi¬ 
tudinal,  lateral  and  normal  accelerations  and  pitch,  roll  and  yaw  rates  as  inputs, 
the  state  variables,  ie.  the  body  axes  velocities  (u,  v  and  w)  and  roll,  pitch  and 
yaw  attitudes  (<p,  9  and  ip)  are  estimated.  From  these  estimates,  the  output  quan¬ 
tities,  velocity,  angle-of-attack,  angle-of-sideslip,  bank  and  pitch  angle  and  altitude 
{y,a,p,<p,0  and  h)  are  calculated  and  compared  with  measured  values  of  the  out¬ 
puts.  The  calibration  constants  are  adjusted  along  with  instrument  bias  parameters 
to  minimise  these  output  errors. 

As  an  example  of  the  procedure  consider  the  force  equation  in  the  Y  body  axis 
direction  for  a  rigid  aircraft  (equation  5. 8, 2(b)  Reference  [11]). 

aerodynamic  contribution  gravity  contribution  inertiat  contribution 

Y  +  mgCosOSintp  =  m(v-t-ru  -pu;)...(l) 

Forces  in  the  X,Y  and  Z  direction  cannot  be  measured  directly  but  use  can  be 
made  of  accelerometers.  Dividing  equation  (1)  by  m 

Y 

- b  gCos9Stn<p  =  (v  -i -  ru  —  pu>)...(2) 
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gives: 


av  +  gCoaQSintp  =  (v  +  ru  -  pu/)...(3) 

The  above  equation  is  a  kinematic  equation  relating  acceleration,  velocity  and  dis¬ 
placement.  Assuming  small  angles  (<  lOdeg.),  a  relationship  between  the  measured 
kinematic  variables  and  angle  of  attack  a  and  angle  of  sideslip  0  can  be  developed: 

a  -  -  yandy  ~  1-(4) 

Then  dividing  both  sides  of  equation  (3)  by  V  gives  : 


*  +  ...(6) 


Substituting  for  a  and  0  gives 

^  CoaBSin<j>  =  0  +  r  -  pa. ..(6) 


Rearranging  equation  (6)  gives  0  in  terms  of  other  kinematic  variables 
0  =  ~  +  ^CoaOSintp  -  r  +  pa. ..(7) 

Now  assuming  that  all  of  the  variables  on  the  right  hand  side  can  be  measured  with 
no  scale  factors  or  bias  errors  0,  can  be  calculated  by  integration  of  equation  (7) 

0  =  J  (^  +  —CoaOSinip  -  r  +  pa)  +  constant. 8) 

A  comparison  can  be  made  between  the  measured  value  of  sideslip  0m  and  the 
value  calculated  from  equation  8  assuming  a  linear  relationship  between  the  local 
Sow  angle  and  the  free  stream  value  of  the  form: 


fim  —  Kff  /^actual  "h 

From  this  comparison  the  scale  factor  Kp  and  bias  bp  can  be  determined. 

For  a  jet  aircraft  with  no  asymmetric  flow,  the  constant  of  integration  or  the 
bias  can  be  expected  to  be  zero.  A  similar  procedure  is  carried  out  with  the  Z  force 
equation  to  determine  an  expression  for  angle  of  attack  a. 

The  estimation  software  used  for  this  purpose  was  developed  for  ARL  under 
a  research  agreement  with  the  University  of  Newcastle.  (References  [12]  and  [13]). 
Application  of  the  method  to  the  measured  time  histories  gave  estimates  for  the  cal¬ 
ibration  constants  for  the  CADS  and  NBTU  angle-of-attack  and  angle-of-sideslip 
measurement  systems.  For  the  angle-of-sideslip  sensor  the  calibration  constant  var¬ 
ied  with  Mach  number.  Values  of  between  1.49  and  1.60  (CADS)  were  used  for 
Phase  1  and  1.06  and  1.20  (NBTU)  for  Phase  2. 
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The  CADS  angle-of-sideslip  sensor  which  is  located  beneath  the  forward  fuselage 
over  estimated  the  true  value  by  50-60%  indicating  strong  cross-flow  in  this  region. 
The  NBTU  gave,  as  expected,  more  accurate  estimates  of  angle-of-sideslip,  over- 
reading  by  only  approximately  10-20%.  While  some  small  variations  occur  in  the 
angle  of  attack  scale  factor,  Ka  with  Mach  number  and  sweep  angle,  these  variations 
are  not  well  defined  within  the  accuracy  of  the  data  and  so  constant  values  of  0.94 
and  1.06  have  been  used  for  Phase  1  and  Phase  2  respectively.  An  investigation 
into  the  effect  of  Ka  variation  showed  that  errors  in  Ka  of  10  %  resulted  in  an 
adjustment  of  typically  5  %  in  the  major  derivatives. 

4.5  .4  priori  Data  from  Model 

A  six  degree  of  freedom  flight  dynamic  model  of  the  F-111C  aircraft,  as  described  in 
Reference  [14),  has  been  developed  at  ARL  and  includes  representation  of  the  flight 
control  system.  A  comprehensive  aerodynamic  data  base  is  used  to  obtain  a  priori  or 
initial  estimates  of  the  aerodynamic  stability  and  control  derivatives.  Configuration 
data  and  initial  conditions  are  defined  for  each  test  point  as  described  in  Appendix 
G. 


4.6  Parameter  Estimation 

A  number  of  techniques  have  been  developed  in  recent  years  for  the  estimation  of 
aerodynamic  derivatives  from  flight  test  measurements.  These  techniques  generally 
use  a  statistical  approach  to  the  process  of  fitting  a  flight  dynamic  model  to  aircraft 
response  time  histories.  The  aerodynamic  stability  and  control  derivatives  are  then 
calculated  from  the  coefficients  of  the  flight  dynamic  model.  For  the  F-111C  data 
analysis,  a  Maximum  Likelihood  technique  was  used.  A  priori  information  for  the 
Maximum  Likelihood  technique  was  obtained  using  the  ARL  six  degree  of  freedom 
flight  dynamic  model  and  is  described  in  Appendix  G.  The  technique  and  associated 
computer  program  are  described  in  References  [15]  and  [16]  and  an  example  MMLE3 
analysis  is  given  Appendix  H. 

Within  the  range  of  flight  conditions  tested  in  the  F-111C  programme,  it  is  as¬ 
sumed  that  the  aircraft  motion  can  be  adequately  represented  by  separate  classical 
linear  flight  dynamic  models  for  longitudinal  and  lateral  motion.  The  longitudi¬ 
nal  and  lateral  flight  dynamic  models,  used  for  this  purpose  are  derived  assuming 
small  disturbance  motion  and  linear  aerodynamic  characteristics.  The  equations 
are  presented  in  Appendix  I  and  J. 

The  maximum  likelihood  technique  uses  a  Newton-Raphson  search  algorithm  to 
iterate  to  a  converged  solution.  The  program  defines  a  cost  function  which  is  the 
weighted  sum  of  the  difference  between  the  model  prediction  and  measured  time 
histories  (see  Figure  7).  Convergence  is  declared  when  the  cost  function  is  less  than 
a  specified  level. 
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The  stabilator  angle  deflection  was  used  as  the  input  for  the  longitudinal  ma¬ 
noeuvres  and  the  rudder,  differential  stabilators  and  spoilers  for  the  lateral  manoeu¬ 
vres.  Note,  spoilers  only  operate  for  cases  with  wing  sweep  less  than  47°.  The  angle 
of  attack,  pitch  rate  and  normal  acceleration  comprise  the  outputs  for  the  longi¬ 
tudinal  model.  The  angle  of  sideslip,  roll  rate,  yaw  rate  and  lateral  acceleration 
comprise  the  outputs  of  the  lateral  model. 

The  derivatives  used  in  the  linear  identification  model  and  their  importance  are 
summarised  in  Figure  3.  The  derivatives  can  be  divided  into  three  categories.  The 
static  derivatives  are  fundamentally  ’stiffness’  parameters  and  the  dynamic  deriva¬ 
tives  are  ’damping’  parameters.  The  analogy  is  drawn  with  a  second  order  mass 
spring  damper  system.  The  control  derivatives  describe  the  forces  and  moments 
due  to  control  surface  deflections.  Figure  9  compares  the  pitch  response  to  an  el¬ 
evator  input  measured  in  flight  with  the  response  calculated  from  the  estimated 
derivatives. 

To  assist  the  identification  procedure,  particularly  for  derivatives  which  make 
only  a  small  contribution  to  the  motion,  a  priori  information  can  be  used  in  the 
identification  procedure.  A  facility  exists  to  constrain  selected  parameters  to  ei¬ 
ther  a  priori  values  or  to  other  model  parameters.  Constraints  were  used  in  the 
longitudinal  model  but  not  in  the  lateral  model. 

For  the  unconstrained  parameters  the  estimation  procedure  calculates  a  measure 
of  the  estimation  accuracy  known  as  the  Cramer-Rao  bound.  The  interpretation  of 
this  quantity  is  given  in  Reference  [15]. 

The  results  plotted  for  each  aerodynamic  derivative  show  an  average  value  for 
this  prediction  error  calculated  from  the  data  points  on  each  plot,  which  can  be  used 
in  conjunction  with  the  observed  repeatability  to  indicate  the  estimation  accuracy. 
To  account  for  the  fact  that  the  signal  noise  is  bandwidth  limited,  the  Cramer-Rao 
bound  is  factored  by  a  multiple  of  10  in  accordance  with  the  procedures  described 
in  Reference  [16].  Two  additional  procedures  were  used  prior  to  the  application  of 
the  maximum  likelihood  technique  to  improve  the  quality  and  consistency  of  the 
measured  time  histories.  These  are  discussed  in  Appendix  H. 

4.7  Curve  Fitting  of  Derivative  Results 

It  is  planned  that  the  data  from  these  flight  tests  will  be  used  to  update  the  aero¬ 
dynamic  data  base  used  in  the  RAAF’s  F-111C  simulator.  To  correct  the  existing 
data-base,  the  new  derivative  information  must  be  presented  as  variations  with 
Mach  number  for  each  sweep  and  altitude.  The  approach  which  has  been  used  is  to 
relocate  the  model  data  curve  to  best  fit  the  flight  derivatives,  therefore  combining 
the  general  information  from  the  flight  tests  and  the  detailed  trends  given  by  the 
model.  The  procedure  used  is  shown  in  Figure  10.  Details  of  this  process  are  pre¬ 
sented  in  Reference  [17],  Details  of  the  notation  used  to  name  analysis  files  is  given 
in  Appendix  K. 
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5  Concluding  Remarks 


This  report  describes  the  steps  involved  in  processing  and  analysing  data  from 
a  flight  test  program  on  an  F-111C  aircraft  to  determine  aerodynamic  stability 
and  control  parameters.  All  stages  in  the  process  are  described  and  details  of 
the  individual  procedures  are  given  in  separate  Appendices.  Computer  programs 
developed  for  the  analysis  are  described  and  running  instructions  are  provided.  The 
programs  and  procedures  can  be  used  for  other  aircraft  test  programmes  with  minor 
alterations. 
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Appendix  A  -  Extracting  Data  From  Flight  Test  Tapes 


The  AFTRAS  program  EXTRACT  is  used  to  extract  flight  test  data  recorded  on 
magnetic  tape.  After  further  processing  by  program  LISTPARM,  a  list  of  parame¬ 
ters  and  associated  values  in  engineering  units  are  outputed  in  ASCII  form.  Each 
channel  has  a  recording  width  of  4000  computer  units  and  calibrations  were  con¬ 
ducted  to  establish  the  relationship  between  computer  units  and  engineering  units. 
A  typical  calibration  is  shown  in  Figure  Al.  Extraction  of  data  at  ARL  is  carried 
out  on  the  EDS  VAX  750  computer  and  the  procedures  involved  are  best  illustrated 
through  the  use  of  an  example  and  reference  to  Figure  A2. 

For  the  purpose  of  this  exercise  the  manoeuvre  selected  occured  on  the  first 
flight  of  the  second  phase  of  tests  and  is  event  72.  From  the  flight  test  report  it  is 
established  that  event  72  commenced  at  approximately  1040.0  seconds  and  finished 
at  about  1109.0  seconds.  From  Table  Al  it  can  be  seen  that  the  data  are  stored  on 
ARL  tape  1153  and  the  data  are  from  the  first  flight  of  a  possible  3  on  this  tape. 

After  logging  onto  the  VAX  750  computer  at  ARL  and  with  tape  1153  mounted 
on  drive  MTB0:,  the  following  commands  are  given  to  assign  devices:  (note  that  all 
user  responses  are  prefixed  by  a  >  symbol) 

>  mount  /foreign/blocksize-16100  mtbO: 

y.MOUNT-1-WRITELOCK.  volume  is  write  locked 
%M0UNT- I -MOUNTED,  mounted  on  _MTB0 : 

>  ass  mtbO:  forOOl 

>  ass  mtbO:  lorOO 2 

>  ass  mtbO:  tape 

Set  default  directory  to  [ae_drobik.aftras.vax. indata]  and  create  the  data  input  file 
P2F1E72.IN.  This  files  contains  the  following  details  of  the  data  to  be  extracted: 

•  The  number  of  channels  to  be  extracted  (33) 

•  Time  interval  (0.01667  or  g^th  of  a  second) 

•  Start  of  extraction  (1040  seconds) 

•  Finish  of  extraction  (1109  seconds) 

•  Channels  and  calibration  versions  to  be  selected  (BI-170  calibration  4  to  BI- 
031  calibration  1) 
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Shown  below  is  file  P2F1E72.IN 


\ 


33 

0. 

01667 

1040.00 

1109.000 

BI 

170 

4 

BI 

171 

4 

BI 

172 

6 

BI 

173 

2 

BI 

174 

3 

BI 

176 

3 

BI 

160 

8 

BI 

201 

7 

BI 

184 

8 

BI 

023 

1 

BI 

026 

2 

BI 

40 

4 

31 

176 

1 

BI 

177 

1 

BI 

214 

4 

BI 

216 

4 

BI 

39 

4 

BI 

41 

1 

CO 

66 

1 

BI 

43 

2 

BI 

44 

3 

BI 

46 

2 

BI 

178 

1 

BI 

63 

2 

BI 

64 

2 

BI 

238 

6 

BI 

190 

8 

BI 

196 

8 

BI 

237 

6 

BI 

188 

4 

BI 

166 

1 

BI 

028 

1 

BI 

031 

1 

Table  A2  lists  the  relevant  channels  in  the  order  to  be  extracted  (1  to  33)  and 
gives  the  calibrations  to  be  used.  The  order  of  the  extracted  channels  is  important 
and  several  points  need  to  be  considered  when  extracting  data: 

•  Channel  7  should  be  selected  according  to  the  altitude  of  the  case  considered. 
For  example  BI- 164-7  would  be  selected  for  a  phase  1  40000  feet  case. 
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•  Channel  8  should  be  selected  according  to  the  Mach  number  of  the  case  con¬ 
sidered  For  example  31-202-8  would  be  selected  for  a  phase  2  Mach  1.2  case. 

•  Channels  10  and  11,  the  angle  of  attack  a  and  angle  of  sideslip  j3  should  be 
selected  according  to  the  phase.  BI-237  and  BI-188  for  phase  1  CADS  a  and 
/?,  or  BI-023  and  BI-025  for  phase  2  NBTU  a  and  j3. 

•  Phase  1  flights  need  only  the  first  28  channels  to  be  extracted. 

•  Phase  2  flights  have  a  total  of  33  channels  with  NBTU  a  and  /?  signals  being 
selected  in  channels  10  and  11. 

•  For  phase  2  flights  channel  31  must  contain  the  crew  module  accelerometer 
BI-31-01  signal.  The  output  is  in  computer  units  and  a  calibration  is  applied 
duririg  the  flight  data  processing  stage.  If  this  channel  is  to  be  used  instead 
of  the  C  of  G  normal  accelerometer  (phase  2  flights  2  and  3)  there  is  no  need 
to  place  this  signal  in  channel  1  because  the  selection  of  normal  accelerometer 
to  be  used  is  carried  out  in  the  FDP  stage. 

Set  default  directory  to  jae.drobik.aftras.vax. extract!  and  run  command  file  EX¬ 
TRACT. COM  as  shown  below: 

>  set  def  [ae_drobik.aftras.vax. extract] 

>  list  extract.com 

3  set  noverify 

$  WRITE  SYS30UTPUT  "  AFTRAS  EXTRACT  PROGRAM  " 

3  WRITE  SYS30UTPUT  "  " 

$  WRITE  SYS30UTPUT  "NB.  Filenames  fully  qualified  eg.  [twp.tsl072lfred.dat" 
$  INQUIRE  q4  "ENTER  THE  MAG  TAPE  DRIVE  e.g.  MTAO:  " 

$  INQUIRE  qi  "ENTER  THE  INput  filename  (.ini)  " 

3  inquire  q2  "ENTER  THE  OUTput  WORKING  filename  (.oul)  " 

3  inquire  q3  "Enter  flight  number  of  the  file  on  the  archive  tape  (1-3) 

II 

3  ru  [ae.drobik. AFTRAS . vax. extract] EXTRACT/PARAMETER3 ( ’Ql ' , ’q2‘ , ’q3 1 . 'q4' ) 

3  WRITE  SYS30UTPUT  "  no  Batch  Job  submitted" 

3  WRITE  SYS30UTPUT  "  -  output  logfile  in  your  root  directory  " 

3  EXIT 

Shown  below  are  the  commands  and  responses  to  commence  the  extract  and  rewind 
the  tape  -a  MTBO: 

>  ru  extract 

F111C  EXTRACT  PROCESSING  PROGRAM 
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Enter  file  code  name,  (eg : P3F2E1Q1)  ::  >  P2F1E72 
Enter  the  tape  sequence  Ho .  ::  >  1 


Inlile- [- . INDATA] P2F1E72 . IN 
Out! lie*  []  P2F1E72 .  OUT 

READING  ARCHIVE  MAG  TAPE 


READ  CALIBRATION  FILE. 

NOW 

CHECKING 

INDATA 

ACCEL  VERT 

II  QM 

BI- 

170- 

4 

ACCEL  LAT 

IIQII 

BI- 

171- 

4 

ACC  LNG  CG 

HQ" 

BI- 

172- 

6 

PITCH  RATE  DEG/SEC 

BI- 

173- 

2 

ROLL  RATE 

DEG/SEC 

BI- 

174- 

3 

YAW  RTE  CG 

DEG/SEC 

BI- 

175- 

3 

HPF  <  6K 

FEET 

BI- 

160- 

8 

MACH  <.7  MACH 

BI- 

201- 

7 

CAS  FINE  T 

KCAS 

BI- 

184- 

8 

NBTU  ALPHA 

DEGREES 

BI- 

23- 

1 

NBTU  BETA 

DEGREES 

BI- 

25- 

2 

STAB  RH  DEGREES 

BI- 

40- 

4 

PTCHACC  CG 

RAD/SEC 

BI- 

176- 

1 

ROLLACC  CG 

RAD/SEC 

BI- 

177- 

1 

ADI  PITCH 

DEGREES 

BI- 

214- 

4 

ADI  BANK 

DEGREES 

BI- 

215- 

4 

STAB  LH 

DEGREES 

BI- 

39- 

4 

RUD  POSN 

DEGREES 

BI- 

41- 

1 

RPEDAL  POS 

INCHES 

BI- 

66- 

1 

SPL  LH  O/B 

DEGREES 

BI- 

43- 

2 

SPL  RH  O/B 

DEGREES 

BI- 

44- 

3 

WNGSWPXDCR 

DEGREES 

BI- 

46- 

2 

YWANG  ACC 

RAD/SEC 

BI- 

178- 

1 

STK  POS  LG 

INCHES 

BI- 

63- 

2 

STK  POS  LT 

INCHES 

BI- 

64- 

2 

AOACADS-VE 

DEGREES 

BI- 

238- 

5 

ALT  CRSCAD 

FEET 

BI- 

190- 

8 

CADS  M  CSE 

MACH 

BI- 

195- 

8 

AOACADS+VE 

DEGREES 

BI- 

237- 

6 

BETA  CADS 

DEGREES 

BI- 

188- 

4 

CM  ACC  VRT 

IIQII 

BI- 

166- 

1 

ACC-X  NBTU 

||  G  " 

51 

28- 

1 

NBTU  A  TMP 

DEGREES 

CBI- 

31- 

1 

FORTRAN  STOP 
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>  SET  MAG/REW  MTBO: 


A  binary  file  of  calibrated  meaaurands  in  engineering  units  has  now  been  cre¬ 
ated.  This  file  has  to  be  converted  to  ASCII  code  and  transfered  to  the  ELXSI 
6400  computer  for  pre  analysis  processing.  The  program  used  for  this  purpose  is 
called  LISTPARM  and  is  invoked  from  the  directory  [ae-drobik.aftras.vax.listparmj. 
Shown  below  are  the  commands  and  responses  to  run  LISTPARM: 

>  RU  LISTPARM 


AFTRAS  LIST  PARAMETERS  PROGRAM 
version  2.1  July  86 
AIRCRAFT  RESEARCH  AND  DEVELOPMENT  UNIT 

PLEASE  ENTER  AFTRAS  OUT  FILENAME 
>  [-. EXTRACT] P2F1E72. OUT 


TYPE  1  IF  PRINTING  BY  EVENTS, 

2  IF  PRINTING  AT  TIME  INCREMENTS, 

3  IF  PRINTING  AVERAGES  BETWEEN  CONSECUTIVE  EVENTS 

4  IF  PRINTING  ‘ABOUT’  EVENTS 

6  IF  PRINTING  ENTIRE  FILE  CONTENTS 

6  OUTPUT  DEVICE  IS  :  SCREEN  ,  SELECT  TO  CHANGE 

7  CHANGE  OUTPUT  MEASURANDS 

8  TO  FINISH 

>  6 

TYPE  1  IF  PRINTING  BY  EVENTS, 

2  IF  PRINTING  AT  TIME  INCREMENTS, 

3  IF  PRINTING  AVERAGES  BETWEEN  CONSECUTIVE  EVENTS 

4  IF  PRINTING  ‘ABOUT’  EVENTS 

6  IF  PRINTING  ENTIRE  FILE  CONTENTS 

8  OUTPUT  DEVICE  IS  :  PRINTER  ,  SELECT  TO  CHANGE 

7  CHANGE  OUTPUT  MEASURANDS 

8  TO  FINISH 

>  6 

FORTRAN  STOP 


Output  from  program  LISTPARM  is  written  to  the  file  LISTPARM.DAT  and 
this  file  is  renamed  to  identify  the  particular  case  extracted  eg.  P2F1E72.DAT. 
This  file  is  copied  to  the  ELXSI  by  logging  onto  the  ELXSI  and  using  the  netcopy 
command. 
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sethost  edsvax 

netcopy  ’edsvax"ae_drobik" : : [ae_drobik.aftras.vax.listparmlP2FlE71.dat 

Flight  data  ia  now  in  ASCII  format  on  the  ELXSI  and  can  be  plotted  using  pro¬ 
grams  TRPLOT  and  TRANS.  Reference  [18]  details  the  TRANS  plotting  program. 
TRPLOT  converts  ASCII  data  to  ELXSI  binary  in  a  format  suitable  for  TRANS 
and  is  run  by: 

>  TRPLOT 

INPUT  DATA  FILENAME  -  >  P2F1E72.DAT 
Read  in  VBLES  etc — beginning  numerical  data 
Finished  with  data--TRFINI  then  end 
Fortran  program  executed  STOP  statement  0 

Now  run  TRANS  to  create  a  device  independent  metafile  plot  file. 

>  TRANS 

[TRANS  version  date  ll-MAR-86] 

I/P  FILENAME  -  plots 

TS  1601  01  A08B-132  006 

I/P  FILE  RECORDED  ON  27-Jul-OO  AT  00:53:64 
INTEGN  INT  =  .0000E+00;  RUN  CPU  TIME  =*  26.84  SEC. 

TIME  FROM  4 . 6800E+02  TO  4.9004E+02  IN  STEPS  OF  1.6670E-02 

>  PLS 

IS  GRAPHICS  OUTPUT  TO  SCREEN  REQUIRED  : 

>  N 

[PLS/O  Output,  lor  this  run,  going  to  DSK:plots.PL6  ] 

STRIP  PLOTS  : 

BLKS 

>  AV 

TO  SPECIFY  NO.  OF  X  UNITS/ INCH.  TYPE  0  FOR  X 
LENGTH  OF  AXES  IN  INCHES;  X.  Y  > 

>  10,6 

ARE  SYMBOLS  REQRD  FOR  PLOTS  : 

>  N 
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LINE  KEY  (0  GIVES  DEFAULT)  => 

>  0 

>  goe 

**  RUNNING  ** 

>  *exi 

A  metafile  PLOTS. PL6  has  now  been  created  and  can  be  viewed  using  a  metafile 
translator  on  any  graphics  device.  Any  channel  can  be  plotted  out  at  the  ARL 
computer  center.  The  file  is  renamed  and  plotted  using  the  following  instructions: 

>  COPY  PLOTS. PL6  PLT.P2F1E71 

>  PLOT . ZT8  PLT.P2F1E72  PICNO-  1  FRAME-634,267 

>  PLOT. ZT8  PLT.P2F1E72  PICNO-  31  FRAME-634,267 


r'* 
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Table  Al:  Data  Storage 


Phase  1  Flights 

|  Flight 

Serials 

ARDU  Tape 

ARL  Tape 

Comments 

No. 

(segment) 

(segment) 

1 

Shake  1 

RDUVX2(1) 

1043(1) 

Shake  down  flight  1 

2 

Shake  2 

RDUVX3(1) 

1044(1) 

Shake  down  flight  2 

3 

1,2 

RXWX2(1) 

1045(1) 

4 

3,4 

RXWX2(2) 

1045(2) 

5 

5,6,7 

RDUVX4(1) 

1046(1) 

6 

10,11 

RDUVX4(2) 

1046(2) 

7 

8,9 

RDUVX4(3) 

1046(3) 

8 

12,13 

RDUVX2(2) 

1043(2) 

9 

14,15 

RDUVX2(3) 

1043(3) 

10 

16 

RDUVX3(2) 

1044(2) 

11 

18,17 

RDUVX3(3) 

1044(3) 

12 

18  repeat 

RDUVX5(1) 

1129(1) 

13 

16  repeat 

RDUVX5(2) 

1129(2) 

14 

17  repeat 

RDUVX5(3) 

1129(3) 

Phase  2  Flights 

Flight 

Serials 

ARDU  Tape 

ARL  Tape 

Comments 

No. 

(segment) 

(segment) 

1 

1  to  10 

JD01(1) 

1153(1) 

2 

11  to  16 

JD01(2) 

1153(2) 

Normal  accelerometer  at 

3 

16,17,18,29,30,31 

JD02(l) 

1067(1) 

CofG  unavailable,  use 

accelerometer  in  crew 

module 

4 

21  to  26,31 

JD02(2) 

1067(2) 

5 

27,28,33,37 

JD04(1) 

1182(1) 

6 

19,20,34,35,36,2  to  7 

JD04(2) 

1182(2) 

7 

27,28,33,37 

JD04(3) 

1182(3) 
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Table  A2:  Instrumentation  Calibrations 


Channel 


Order 

Channel  Calibration  Number 

Phase  1 

Phase  2 

1 

444444444444 

4444444 

2 

444444444444 

4444444 

3 

445555555555 

5555555 

4 

222222222222 

2222222 

5 

333333333333 

3333333 

6 

333333333333 

3333333 

acceleration  vert. 

170 

acceleration  lat. 

171 

acceleration  long. 

172 

Pitch  rate 

173 

Roll  rate 

174 

Yaw  rate 

175 

Height  <  5  k 

160 

Mach  <  0.7 

201 

CAS  fine 

184 

AoA  CADS  +ve 

237 

Beta  CADS 

188 

Right  stab 

040 

Pitch  acceleration 

176 

Roll  acceleration 

177 

ADI  pitch 

214 

ADI  bank 

215 

Left  stab 

039 

Rudder  position 

041 

Rudder  pedal  posn. 

055 

Spoiler  LH  O/B 

043 

Spoiler  RH  O/B 

044 

Wing  sweep 

045 

Yaw  acceleration 

178 

Long,  stick  posn. 

053 

Lat.  stick  posn. 

054 

AoA  CADS  -ve 

238 

Alt  coarse 

190 

CADS  Mach  coarse 

195 

666666666666 

666666666666 


111111111111 

111111111111 


111111111111 

111111111111 

222222222222 

333333333333 

222222222222 

111111111111 


8  8  8  8  8  8  8 


7777777 

8888888 


1111111 

1111111 


1111111 

1111111 


1111111 
2  2  2  2  2  2  2 
2  2  2  2  2  2  2 
5  5  5  5  5  5  5 
8888888 
8888888 


Additional  Channels  for  Phase  2  Flights 


0000000 


Available  Height  and  Mach  Range  Channels 


Height  <  5k 
5k  <  Height  <  15k 
15k  <  Height  <  25k 
25k  <  Height  <  35k 
35k  <  Height  <  45k 
45k  <  Height  <  55k 
Mach  <  0.7 
0.69<  Mach  <  1.8 


777777777777 

777777777777 

777777777777 

777777777777 

666666666666 

666666666666 

666666666666 

777777777777 


8  8  8  8  8  8  8 
8888888 
8888888 
8888888 
7777777 
7  7  7  7  7  7  7 
7777777 
8888888 
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Figure  A2:  Extracting  Flight  teat  data 
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Appendix  B  Mass  Characteristics  Program  CGCALCP10R3 

Computer  program  CGCALC.P10R3  La  written  in  the  computer  language  FOR¬ 
TRAN  77  and  uses  manufacturers  and  experimental  data  to  calculate  the  aircraft’s 
mass  characteristics.  Figure  B1  gives  an  outline  of  the  program.  At  run  time  the 
user  has  a  number  of  options  available  including: 

•  Phase  1  or  phase  2  flight 

•  CADS  or  NBTU  a  and  /?  measurement  to  be  used  for  analysis.  This  infor¬ 
mation  is  required  to  allow  calculation  of  the  distance  between  the  location  of 
the  transducers  and  the  C  of  G. 

•  Source  of  normal  accelerometer  measurement. 

•  A  normal  accelerometer  signal  was  available  from  3  sources,  C  of  G,  crew 
module  or  from  the  nose  boom  in  phase  2  flights  and  like  the  a  and  f)  signals, 
its  position  relative  to  the  C  of  G  has  to  be  calculated. 

•  Slat  retracted  or  extended 

•  Flap  extension.  Note  that  if  flaps  are  extended,  the  slats  will  also  be  extended. 

•  Method  of  fuel  contents  calculation:  Either  total  fuel  contents  along  with  the 
automatic  fuel  schedule,  or  the  indicated  fuel  contents  of  the  forward  and 
aft  tanks  corrected  with  fuel  calibration  data,  to  calculate  actual  fuel  tank 
contents.  Table  Bl  shows  the  automatic  fuel  schedule  followed  if  this  option 
is  used. 

A  typical  run  of  program  CGCALC-P10R3  is  shown  below  and  user  inputs  are 
indicated  as  > 

>  CGCALC.P10R3 

WING  SWEEP  ANGLE  IS 

>  16 

PHASE  1  OR  2  FLIGHT  [1  OR  2] 

>  2 

SPECIFY  IF  CADS  OR  NBTU  * ALPHA*  MEASUREMENT  TO  BE  USED 
[RESPONSE  IS  EITHER  CADS  OR  NBTU] 

>  NBTU 

SPECIFY  IF  CADS  OR  NBTU  *BETA*  MEASUREMENT  TO  BE  USED 
[RESPONSE  IS  EITHER  CADS  OR  NBTU] 

>  NBTU 

SPECIFY  WHERE  THE  NORMAL  ACCELEROMETER  IS  LOCATED  AT 
THE  CofG  ,  ON  THE  NBTU  OR  IN  THE  CREW  MODULE 
[RESPONSE  IS  CG  ,  NBTU  OR  CREW] 
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>  CREW 

ARE  SLATS  EXTENDED  [RESPNSE  IS  YES  or  NO] 

>  YES 

FLAP  DEFLECTION 

>  26.0 

DO  YOU  WANT  TO  SPECIFY  INDICATED  FUEL  TANK  CONTENTS 
FOR  AIRCRAFT  A8-132  [Y  OR  N] 

>  Y 

FORWARD  FUEL  TANK  READING  IS  [LBS] 

>  14300 

AFT  FUEL  TANK  READING  IS  [LBS] 

>  6000 


AIRCRAFT  CONFIGURATION  WITH  INDICATED  FUEL 


WEIGHT 

DELTA  X 

DELTA  Y 

DELTA  Z 

BASIC  AIRCRAFT 

40073.6 

520.1 

.0 

180.2 

SWEPT  WING 

7174.4 

623.7 

40.0 

176.2 

PILOT 

216.0 

249.7 

-13.0 

184.2 

NAVIGATOR 

216.0 

240.7 

13.0 

184.2 

OIL  (UNUSEABLE) 

09.0 

686.9 

.0 

164.6 

OIL  (USEABLE) 

12.0 

686.9 

.0 

168.3 

WATER  (AIR  CON) 

41.7 

470.0 

.0 

166.0 

PYROTECHNICS 

24.8 

801.1 

.0 

192.0 

OXYGEN  CONVERTER 

26.0 

162.8 

.0 

167.0 

CONVERTER 

17.0 

162.8 

.0 

167.0 

CREW  MODUALE  BALLAST 

27.0 

268.0 

.0 

170.0 

WINDSHEILD  WASH  FLUID 

8.1 

264.0 

.0 

143.0 

EMG  OXY  AND  COMP  AIR 

6.8 

282.9 

.0 

207.6 

MISSION  DATA  PROV. 

3.0 

267.9 

.0 

179.0 

CHAFF 

37.0 

801.1 

.0 

192.0 

GUN  (EMPTY) 

1061.8 

386.6 

.0 

162.0 

FORWARD  TANKS 

14449.3 

400.6 

.0 

191.6 

AFT  TANKS 

6933.6 

634.8 

.0 

179.9 

WING  TANKS  TOTAL 

.0 

616.9 

.0 

194.8 

WEAPON  BAY  TANK 

.0 

392.0 

.0 

146.4 

FUEL  IN  LINE 

184.0 

616.5 

.0 

184.0 

UNUSEABLE  FUEL 

292.0 

494.5 

.0 

179.3 

SPARE1 

.0 

.0 

.0 

.0 

SPARE2 

.0 

.0 

.0 

.0 

SPARE3 

.0 

.0 

.0 

.0 

FWD  TANKS  CONTAIN  8616.72363  POUNDS  OF  FUEL  MORE  THAN  AFT 
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FORWARD  FUEL  GAUGE  READING3  14300.0  POUNDS 
ACTUAL  FORWARD  FUEL  -  14449.291  POUNDS 

AFT  FUEL  GAUGE  READING  3  6000.0  POUNDS 
ACTUAL  AFT  FUEL  3  6933.66738  POUNDS 

FUEL  USED  3  14199.1406  POUNDS 

HQRZ.  POSITION  OF  CG  3  24.4684624  PERCENT  OF  MAC 

HORZ.  POSITION  OF  CG  3  607.008697  FUSESTAX 

VERTICAL  POSITION  OF  CG  3  3.83693628  PERCENT  OF  MAC 

VERTICAL  POSITION  OF  CG  3  181.661803  WATERLNZ 

TOTAL  USEABLE  FUEL  WEIGHT  >20382. 8693  POUNDS 
TOTAL  WEIGHT  OF  AIRCRAFT  >69900.8693  POUNDS 


FLAP  DEFLECTION  =26.0 

V******  INERTIA  DATA  ****** 


IXX 

>  73343.6718 

SLUGS  FT-SQ 

IYY 

>  344866.781 

SLUGS  FT-SQ 

IZZ 

3  411017.062 

SLUGS  FT-SQ 

IXZ 

3  3976.7871 

SLUGS  FT-SQ 

-IXZ/IXX  -  -6.42212724E-02 


-IXZ/IZZ  -  -9.67647949E-03 

**********  INSTRUMENT  OFFSET  DATA  ********** 


XALF 

3 

46.6107 

YALF 

3 

.0000 

ZALF  3 

2.6636 

XB 

3 

44.7032 

YB 

3 

.0000 

ZB  3 

3.0618 

XAX 

3 

2.5241 

YAX 

= 

-2.4308 

ZAX  3 

1 . 3877 

XAY 

3 

2.4218 

YAY 

3 

-2.4583 

ZAY  3 

1.3877 

XAN 

m 

20.8674 

YAN 

3 

.0000 

ZAN  3 

1.6610 
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Table  Bl:  Automatic  Fuel  Schedule 


TANK 

Fuel  Used  In 
Segment  (lbs) 

Total  fuel  used 
(lbs) 

Fuel  to  Burn 
(lbs) 

Full  Fuel  Load 

0 

0 

34582 

Primary  Weapon 

Bay  Fuel  Tank 

1852 

1852 

32730 

Wing  tanks 

5060 

2760 

27670 

Forward  1  tank 

712 

7624 

26958 

Forward  1  tank  it 

Aft  2  tank  until 

Aft  2  tank  empty 

4868 

12492 

22090 

Forward  1  tank  & 

Aft  1  tank  until 
Forward  1  tank  empty 

9142 

21634 

12948 

Forward  2  tank  it 

Aft  1  tank  until 

Aft  1  tank  empty 

4564 

26198 

8384 

Forward  2  tank 
until  empty 

5743 

31941 

Reserve 

2457 

34398 

184 

Fuel  in  Lines 

184 

34582 

0 
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User  input  Induces 


1)  Wing  sweep 

2)  alpha  and  beta  CADS  or  NBTU 

3)  Normal  accelerometer  location 

4)  Slat  and  flap  extensions 

5)  Indicated  fuel  and  use  calibration  or 
total  fuel  and  use  fuel  scheduie 


Baseline  aircraft  data 


Data  base  containing  l  of  G,  Inertia, 

flap,  slat  extension  information 

V _ / 


Program  CGCALC _ P10R3 


Calculate  individual  fuel  tank  contents  if  using  fuel  schedule  or  i 
apply  calibration  to  indicated  fuel  tank  readings  j 


Calculate  C  of  G  of  empty  wing 


Calculate  C  of  G  of  following  fuel  tanks:  Wing 

Primary  weapon  bay 

Forward 

Aft 


Calculate  C  of  G  change  if  slats  and  flaps  are  extended 


Adu  component  contnbutions  to  basic  aircraft  weight  and 
C  of  G  and  calculate  overall  AUW  and  C  of  G 


Calculate  instrumentation  reference  data  ] 


|  Calculate  moments  and  products  of  inertia  for  various  components 
I  about  their  own  C  of  G  and  then  calculate  overall  values 


OUT01  JT: 

Files  required  for  use  in  Flight  Data  Processing 
and  parameter  estimation 


Figure  Bl:  Mass  Characteristics  Calculation 
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Appendix  C  -  Flight  Data  Processing  Program  FDP 

The  flight  data  processing  program  FDP  reads  in  flight  data  obtained  from  the 
AFTRAS  program  LISTPARM  and  processes  it  into  files  formatted  for  use  in  a 
variety  of  analysis  software  programs.  Using  the  manoeuvre  given  in  Appendix  B, 
an  example  of  running  FDP  is  shown  below.  A  flow  diagram  showing  the  program 
steps  is  given  in  Figure  Cl.  (note  that  all  user  responses  are  prefixed  by  a  >): 

>  FD 

****  WELCOME  TC  F111C  FLIGHT  DATA  PROCESSING  PROGRAM  -  FDP  **** 
*******  RUNNING  PROGRAM  FDP  TO  PROCESS  FLIGHT  DATA  ********* 

Enter  listparm  input  *ile  name  in  lull 

>  P2F1E72.DAT 

Using  the  instrument  position  information  contained  in  the  file  MMLE3.PARA, 
which  is  an  output  from  program  CGCALC_P10R3,  the  instrument  options  are 
c  is  layed  and  can  be  used  to  check  that  all  required  options  are  correct.  If  the  case 
to  ok  processed  is  from  a  phase  2  flight,  the  phase  2  pressure  error  correction  will 
be  applied.  Also  printed  out  are  the  start  and  finishing  times  of  the  extracted  data 
and  the  data  channels  and  their  respective  calibrations.  A  check  should  be  made 
at  this  stage  to  ensure  that  the  required  channels  have  been  extr^ted. 

NBTU  Alpha  instrument  is  being  used 

NBTU  Beta  instrument  is  being  used 

Crew  module  Normal  Accel  Instrument  is  being  used 

Phase  2  PE  Corrections  Applied 

INFORMATION  ABOUT  DATASET 


NO.  OF  INC.  START  TIME  END  TIME 

MEASURANDS 

33  1 . 668Q00O8E-O2  1040.0  1100.0 

DATASET  DESCRIPTION 

TRIAL  ID  TRIAL  PHASE  A/C  ID  FLIGHT  NO. 

TS  1601  02  A08B-132  001 

NO .  SHORT  SHORT  MEASURAND 
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TITLE 

UNITS  DESCRIPTION 

1 

ACCEL  VERT 

"G" 

BI-170-  4 

2 

ACCEL  LAT 

HQ" 

BI-171-  4 

3 

ACC  LNG  CG 

"G" 

BI-172-  6 

4 

PITCH  RATE 

DEG/SEC 

BI-173*  2 

6 

ROLL  RATE 

DEG/SEC 

BI- 174-  3 

6 

YAW  RTE  CG 

DEG/SEC 

BI-176-  3 

7 

HPF  <  6K 

FEET 

BI-180-  8 

8 

MACH  <.7 

MACH 

BI-201-  7 

9 

CAS  FINE  T 

KCAS 

BI- 184-  8 

10 

NBTU  ALPHA 

DEGREES 

BI-  23-  1 

11 

NBTU  BETA 

DEGREES 

BI-  26-  2 

12 

STAB  RH 

DEGREES 

BI-  40-  4 

13 

PTCHACC  CG 

RAD/SEC 

BI-176-  1 

14 

ROLLACC  CG 

RAD/SEC 

BI- 177-  1 

16 

ADI  PITCH 

DEGREES 

BI-214-  4 

16 

ADI  BANK 

DEGREES 

BI-215-  4 

17 

STAB  LH 

DEGREES 

BI-  39-  4 

18 

ROD  POSN 

DEGREES 

BI-  41-  1 

19 

RPEDAL  POS 

INCHES 

BI-  66-  1 

20 

SPL  LH  0/B 

DEGREES 

BI-  43-  2 

21 

SPL  RH  0/B 

DEGREES 

BI-  44-  3 

22 

WNGSWPXDCR 

DEGREES 

BI-  46-  2 

23 

YWANG  ACC 

RAD/SEC 

BI-178-  1 

24 

STK  POS  LG 

INCHES 

BI-  63-  2 

26 

STK  POS  LT 

INCHES 

BI-  64-  2 

26 

AOACADS-VE 

DEGREES 

BI-238-  6 

27 

ALT  CRSCAD 

FEET 

BI- 190-  8 

28 

CADS  M  CSE 

MACH 

BI-195-  8 

29 

AOACADS+VE 

DEGREES 

BI-237-  6 

30 

BETA  CADS 

DEGREES 

BI-188-  4 

31 

CM  ACC  VRT 

11 G" 

BI-166-  1 

32 

ACC-X  NBTU 

"G» 

BI-  28-  1 

33 

NBTU  A  TMP 

DEGREES  C 

BI-  31-  1 

The  option  is  available  at  this  stage  to  select  different  portions  of  the  extracted 
time  histories  for  analysis. 

Do  you  wish  to  delete  data  from  start  of  time  history  ? 

>  Y 

Enter  Aftras  tape  time  to  begin  time  history 
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>  1089 


Do  you  wish  to  delete  data  irom  end  of  time  history  ? 

>  Y 

Enter  Altras  tape  time  for  data  cut-oli 

>  1080 

Type  of  phase  lags  to  be  used 


Do  you  want  standard  Channel  lags  ? 

For  standard  Phase  1  lags  type  . [PI] 

For  standard  Phase  2  lags  type  . [P3] 

For  standard  Phase  2  lags  with  the  crew 

module  norm  accelerometer  type  . [CM] 

For  non  standard  lags  type  . [N] 


>  CM 

Standard  Phase  2  Crew  Module  lags  applied 
**  Flight  Data  now  being  processed  ** 


The  processing  of  the  raw  flight  data  now  commences  and  several  files  are  created 
for  use  in  the  analysis  as  shown  below: 

•  FPR.LON  and  FPR.LAT.  These  files  provide  the  input  data  for  the  flight 
path  reconstruction  programs  ARLFPRLON2  and  ARLFPRLAT  which  are 
used  to  calculate  the  flow  vane  calibration  constants.  Note:  the  data  for  these 
programs  are  required  in  the  following  units.  Lengths  in  meters,  angles  in 
degrees,  angular  rates  in  radians  per  second,  velocity  (VTAS)  in  meters  per 
second  and  accelerations  at  the  centre  of  gravity  in  meters  per  second  per 
second.  Because  the  instrumentation  package  was  not  located  at  the  centre 
of  gravity,  corrections  have  to  be  made  to  the  measured  signals  to  give  the  C 
of  G  values.  Using  Figure  C2  and  the  equations  given  below  the  C  of  G  offset 
corrections  are  made. 


A  zq  Ayr  A  xq1  Air2 

- I - 1 - -j - 

9  9  9  9 
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av«,  =  av 


A  xr  A  zp  Ayr2  Ayp2 

- 1 - H - i - 

9  9  9  9 

A  xq  Ay  p  Azg2  Azp2 

9  T  9  9  9 


•  MMLE3.DAT.  This  is  the  raw  data  file  required  for  use  in  the  parameter 
identification  program  MMLE3.  The  data  for  this  program  are  required  in 
the  following  units: 

angles  in  degrees,  angular  rates  in  radians  per  second  velocity  (VTAS)  and 
accelerations  in  ’g’  measured  at  that  instrumentation  pack. 


•  F-lllC.dat.  This  file  can  be  used  to  provide  control  input  time  histories 
into  a  six  degree  of  freedom  model  of  the  Flll-C  enabling  comparisons  to 
be  made  between  the  calculated  response  of  the  aircraft  during  the  flight  test 
manoeuvre  and  the  actual  measured  response. 


Have  a  glorious  day  in  Melbourne  I 


************  flight  data  processing  complete  **************** 

The  flight  data  processing  has  been  completed  and  analysis  can  commence. 
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Figure  Cl:  Flight  Data  Processing 


37 


Appendix  D  -  Pressure  Error  Correction 

For  the  manoeuvres  flown  with  the  NBTU  fitted  a  correction  is  required  to  the  data 
calculated  by  the  Central  Air  Data  System  due  to  the  change  in  the  static  pressure 
source.  This  correction  has  to  include  the  differences  between  the  static  pressures 
at  the  two  locations  and  also  the  in-built  corrections  which  are  used  by  the  CADS 
system  for  data  from  the  aircraft  system.  A  calibrated  Mirage  chase  aircraft  was 
used  to  establish  the  magnitude  of  the  error  of  the  static  pressure  source  for  aircraft 
A8-132  with  the  NBTU  fitted.  The  static  and  impact  pressures  can  be  calculated 
from  the  following  equations  taken  from  Reference  [19].  The  static  air  pressure  p,p 
for  the  pacer  Mirage  aircraft  at  a  level  hpaeer  is  given  by: 

(  To  —  Lhpaccr  \  A 

=  K - 5b - ) 

The  ambient  pressure  pa  can  be  then  calculated  using  the  known  pressure  error 
correction  from  the  pacer  aircraft. 

For  the  F-111C  aircraft  with  an  indicated  altitude  of  h,  the  static  pressure  p,  is 
calculated  from: 

/T0- Lhi\T?n 
P.  Po[  Tq  ) 

The  differential  pressure  A p  is 


Ap  =  p,  -  Pa 


From  the  F-lllC’s  indicated  airspeed,  VIAS  the  indicated  impact  pressure  qci  is 
calculated 


—  P  o 


'r  (7-1)  /viAsy 

7-1 

-  1 

2  V  a0  / 

The  chase  and  F-111C  aircraft  were  flown  in  close  formation  at  various  values 
of  VIAS  and  the  indicated  altitude  from  the  two  aircraft  was  recorded. 

This  data  was  analysed  by  ARDU  using  the  known  pressure  correction  of  the 
chase  aircraft  to  give  the  static  pressure  corrections  shown  in  Table  Dl.  Measure¬ 
ments  were  not  made  at  supersonic  speeds.  Data  from  the  subsonic  tests  were 
extrapolated  to  supersonic  speeds  using  information  from  tests  on  similar  installa¬ 
tions. 

The  corrections  to  velocity  and  altitude  are  calculated  using  the  values  of  AP/qci 
given  in  Table  Dl  by  the  following  method: 

For  the  given  indicated  pressure  height  hi  calculate  the  static  pressure  p,  using 
the  following  equations: 

For  H{  <  36089.0  feet 
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For  Hi  >  36089.0  feet 

Pt  =  pie*7v 

where  T\  and  pi  are  the  temperature  and  static  pressure  at  36089  feet 
Given  the  indicated  airspeed  VLAS,  the  indicated  impact  pressure  is  calculated 
using  the  following  equation  : 


9ci  =  Po 


l  (7-1)  /V/A5\2] 

^  1 

2  \  ao  / 

Using  qd  and  p,  calculated  above,  the  indicated  Mach  number  M{  is  calculated: 


M,= 


7-1  Vp,  ) 


From  Table  Dl  for  any  given  value  of  Mi  interpolate  for  a  value  of  A p/qci  from 
which: 

A p  =  A p/qci  x  qd 

can  be  calculated.  The  actual  impact  pressure  is  calculated  to  be  : 


9c  =  qd  +  Ap 


Using  this  calibrated  impact  pressure  calculate  the  calibrated  airspeed  VCAS 


VC  AS  =  a0 


0.5 


The  ambient  calibrated  pressure  is  : 


p9  =  p,  -  A p 


And  the  calibrated  pressure  altitude  is  : 

For  Hi  <  36089.0  feet 

h ” - r 

For  Hi  >  36089.0  feet 

hpe  =  ~ log ,  (— )  +  36089 
9  \Pi/ 

Where  T\  and  pi  are  the  temperature  and  static  pressure  at  36089  feet 
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Table  Dl:  Pressure  Error  Correction  Data 


Mach  No. 

Ait 

V,' 

?c, 

A  P 

A  P/qci 

1000 

325 

379.71 

16.127 

0.0425 

1000 

390 

561.26 

28.956 

0.0516 

BE 

1000 

456 

791.64 

45.215 

0.0571 

1000 

520 

1066.0 

68.611 

0.0644 

m 

1000 

588 

1420.0 

112.000 

0.0787 

BM 

1000 

620 

1613.0 

145.000 

0.0902 

T*- 

\ 


Appendix  E  -  Scale  Altitude  Correction 


Scale  altitude  error  is  the  term  used  to  describe  the  difference  between  calibrated 
airspeed  VCAS,  which  is  the  indicated  airspeed  VIAS  corrected  for  any  pressure  er¬ 
rors,  and  instrumentation  errors,  and  the  equivalent  airspeed  VEAS.  The  difference 
between  VCAS  and  VEAS  occurs  because  airspeed  measuring  instrument  calibra¬ 
tions  are  simplified  by  assuming  a  constant  value  of  unity  for  the  pressure  ratio 
6  =  p/po  in  part  of  the  calibration  equation.  This  error  can  be  neglected  at  low 
speeds  and  low  altitudes  but  is  significant  at  high  Mach  numbers  and  high  altitudes. 
Bernoulli’s  equation  for  the  total  pressure  in  compressible  flow  can  be  expressed  as 
a  function  of  Mach  number  in  the  absence  of  any  shock  wave  as: 


Pt  =  P 


i  + 

2 


i-i 


...{El) 


which  is  the  pitot  pressure  pp  sensed  by  a  pitot  tube  in  subsonic  flow.  For  supersonic 
flow,  the  pitot  pressure  is  modified  by  the  presence  of  a  normal  shock  at  the  mouth 
of  the  pitot  tube.  The  ratio  of  total  pressure  across  a  shock  is  given  by  : 


Pti  = 
Pti 


1  + 


27 


(M2  -  1) 


(1  +  1) 

Therefore  for  M  <  1.0,  pp  =  pt 

(l_-  A! 


l“1 


,..{E2) 


P 


I  +  -L~— M 
2 


i-i 


...(E3) 


and  for  M  >  1.0,  pp  ^  pt 


Pp  _  P*2_  x  Pfi_  _ 
P  Pi,  P 


1  + 


27 


(7  +  1) 


(M2  -  1) 


T+t  r 


(x+Hm2 

2 


i-i 


,(E4) 


The  impact  pressure,  which  is  the  difference  between  pitot  and  stagnation  pressure 
is  given  by  the  relationship  : 


9c  =PP  ~  p.  .{E5) 

Calibrated  airspeed  VCAS  is  obtained  from  the  indicated  airspeed,  corrected  for 
any  pressure  measurement  errors  that  have  been  introduced  by  the  air  data  system 
and  for  any  mechanical  errors  in  the  instrument.  VCAS  is  the  airspeed  that  would 
be  required  at  sea-ievel  in  a  standard  atmosphere  to  give  the  same  impact  pressure 
as  that  sensed  at  the  particular  height  and  true  airspeed  considered.  For  a  sea  level 
ISA  atmosphere: 

When  VCAS/ao  <  1.0 


9c  _  Pp  ~  P 
Po  Po 


+ 


(7~1) 

2 


VCAS\ 
a  o  ) 


l...(£6) 
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When  VCAS/ao  >  1.0 


?C  _  Pp  ~  P  _ 


PO  PO 


=  1  + 


(7+1) 


To  find  the  relationship  between  calibrated  airspeed  and  Mach  number  we  can 
equate  the  above  expressions  using  the  substitution 


Note:  Three  cases  have  to  be  considered:  For  VCAS/ao  <  1-0  and  M  <  1.0 
equating  E6  and  E3 

'f  +h-ii/V£i£y|*_l  fr 

2  V  a0  /  [L  2 

For  VCAS/ao  <1.0  and  M  >  1.0  equating  E6  and  E4 

■[  b-o  /vc^yi* _  1  , fr,  +  *.  (M, . 1}|^ _ ,1 

2  V  ao  )  [L  (7+1)  J  (2 


For  VCAS/ao  >1.0  and  M  >  1.0  equating  E7  and  E4 


(7  +  1)  V  V  ao 


21  TT^TT 


2q  ( fVCAS\2  ,\]^\(1  +  l)  fVCAS\2]^  , 

;  (7+1)  VV  ao  )  VJ  2  V  ao  )  . 

Because  6  is  a  function  of  pressure  height,  ie  p  is  the  ambient  pressure  at  a  given 
height,  the  above  equations  show  that  there  is  a  relationship  between  VCAS,  M  and 
pressure  height.  Given 

VEAS  =  VTAS(-) 2 

Po 


VEAS  =  Ma 


(£>* 
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=  Mao6 3 

Substituting  for  VEAS  =  aoMS »  and  7  =  1.4  the  relationships  between  VCAS, 
VEAS  and  pressure  height  can  be  developed  for  the  three  cases. 

For  VCAS/ao  <  1.0  and  VEAS/{a0 7)0  6  <  1.0  then 


For  VCAS/ao  <  1.0  <  V C  AS / (ao'r)0  s  then 


1+0.2 


/ VCAS\ 2 

V  ao 


)“-*] 

II 

o> 

l 

12**6 


hi  (vEAsy 


[  (7  6  (kJeas)  ) 


2.5 


-  1 


For  VCAS/ao  >  1.0  and  VCAS/(a0 7)0  6  >  1.0  then 


1  28  x  52  6  (V£M)2 

(7  -  (r*B)’)“  * 


=  6 


1  2»xSai>  (VEAS\: 
6  \  00  ) 

(7_5(Ffe)2) 


-  1 


From  the  above  expressions  the  difference  between  calibrated  airspeed  VCAS  and 
the  equivalent  air  speed  VEAS  can  be  calculated.  At  sea  level  6  =  1  so  the  values 
of  VCAS  and  VEAS  are  identical.  At  higher  altitudes  where  6  ^  1.0  the  difference 
between  the  two  increases.  This  difference  has  been  called  scale  altitude  error. 
Figure  El  shows  the  value  of  scale  altitude  error  as  a  function  of  pressure  height 
and  Mach  number.  For  airspeed  measurement  purposes  it  is  necessary  to  make  an 
adjusment  to  VCAS  especially  at  high  altitudes  and  Mach  numbers. 
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Figure  EL:  Scale  Altitude  Correction  to  VCAS 


Appendix  F  -  Flight  Path  Reconstruction  Programs 

The  programs  used  are  described  in  detail  in  References  [12]  and  [  13]  and  in  this 
section  an  example  will  be  used  to  illustrate  how  they  are  run  on  the  ELXSI  com¬ 
puter.  The  first  step  is  to  set  up  the  data  input  table  ’DATIN.DAT’  and  this  is 
done  using  program  FPRSORT.  As  with  previous  examples  user  input  is  prefixed 
by  a  >  symbol. 

>  fprsort 

NAME  OF  INPUT  FILE 

>  datin.long 

LONGITUDINAL  CASE 
TABLE  1  :  STATE  VARIABLES 

1  U  YES 

3  W  YES 

6  THE  YES 

7  H  YES 

9  BAY  NO 

11  BP  NO 

13  BR  NO 

16  LAY  NO 

17  LP  NO 

19  LR  NO 

21  BBET  NO 

23  BPHI  NO 

26  BPSI  NO 

27  LVEL  NO 

29  LA LP  YES 

31  LTHE  NO 

33  LH  NO 

TABLE  2  :  INPUTS 

1  Ax  YES 

3  Az  YES 

6  q  YES 

TABLE  3  :  OUTPUTS 


1 

V 

YES 

2 

BET 

NO 

3 

ALP 

YES 

4 

PHI 

NO 

6 

THE 

YES 

6 

PSI 

NO 

7 

H 

YES 

8 

2  Ay  NO 

4  P  NO 

6  R  NO 


2 

V 

NO 

4 

PHI 

NO 

6 

PSI 

NO 

8 

BAX 

YES 

10 

BAZ 

YES 

12 

Bq 

YES 

14 

LAX 

NO 

16 

LAZ 

NO 

18 

Lq 

NO 

20 

BVEL 

YES 

22 

BALP 

YES 

24 

BTHE 

YES 

26 

BH 

YES 

28 

LBET 

NO 

30 

LPHI 

NO 

32 

LPSI 

NO 

34 
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TABLE  #  TO  BE  EDITED? 

(0  WILL  TERMINATE  CHANGES, -1  WILL  GIVE  TABLE  CONTENTS): 

The  state  variables,  inputs  and  outputs  to  be  used  for  analysis  are  displayed 
and  can  be  changed.  7  „r  this  case  (a  longitudinal  pitch  up  manoeuvre),  the  option 
for  changing  the  states  variable  inputs  and  outputs  are  not  exercised. 

>  0 

The  start  and  finishing  times  are  displayed  along  with  the  position  (in  meters)  of 
the  angle  of  attack  a  and  angle  of  sideslip  measurement  devices.  In  the  listing 
below  the  finishing  time  TEND  is  changed  to  6.49  seconds  and  the  distance  to  the 
a  measuring  device  XALP  is  increased  to  13.862  meters. 


1 

TO 

.  OOOOOE+OO 

2 

TEND 

.40000E+02 

3 

DT 

. 60000E+02 

4 

XALP 

. 46000E+01 

6 

YALP 

. 00000E+00 

6 

XBET 

.OOOOOE+OO 

7 

ZBET 

.OOOOOE+OO 

NO.  AND  VALUE  TO  CHANGE  (EG.  3  .01)  OR  0  0  TO  TERMINATE 
2  6 . 49 


1 

TO 

. OOOOOE+OO 

2 

TEND 

.64000E+01 

3 

DT 

. 80000Z+02 

4 

XALP 

. 46000E+01 

6 

YALP 

.OOOOOE+OO 

6 

XBET 

. OOOOOE+OO 

7 

ZBET 

.OOOOOE+OO 

NO.  AND  VALUE  TO  CHANGE  (EG.  3  .01)  OR  0  0  TO  TERMINATE 
4  13.862 


1 

TO 

.OOOOOE+OO 

2 

TEND 

.64000E+01 

3 

DT 

. 60000E+02 

4 

XALP 

. 13862E+01 

6 

YALP 

.OOOOOE+OO 

6 

XBET 

.OOOOOE+OO 

7 

ZBET 

.OOOOOE+OO 
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5 
0 

154 
1 
1 
1 
1 

THE  LAST  FOUR  PARAMETERS  ARE  TRUE  OR  FALSE  SWITCHES 
0  -  FALSE  1  -  TRUE 
0  0 

For  longitudinal  manoeuvres  the  program  ARLFPR2LON  is  used  and  for  lateral 
cases  program  ARLFPRLAT  is  used.  For  the  program  to  commence  running  type: 

>  arlfpr21on 

Shown  below  is  part  of  the  files  DATOUT.DAT  which  summarises  the  results  ob¬ 
tained: 


NO.  AND  VALUE  TO  CHANGE  (EG.  3  .01)  OR  0  0  TO  TERMINATE 
0  0 

1  NPAS  No  of  passes  through  data 

2  ITER  No  of  local  Iterations 

3  NDEI  No  of  samples  before  averaging 

4  IFIL  filtering  of  inputs  write  to  FILINP.DAT 

5  ICYL  samples  between  writes  to  the  output  files 

0  IRES  residuals  k  S.D.  written  to  RES. DAT  k  SDEV.DAT 

7  IBIAS  inst.  param.  written  to  PARAM.DAT  on  last  pass 

NO.  AND  VALUE  TO  CHANGE  (EG.  3  3800)  OR  0  0  TO  TERMINATE 


**************************************************************************** 
FLIGHT  PATH  RECONSTRUCTION:  PROGRAM  ARLFPR  RUN  1 
**************************************************************************** 

**  PARAMETERS  OF  THE  FILTERING  PROCEDURE  ** 

START  TIME  «  .00000E+00  SECS  **  FINAL  TIME  -  .64000E+01  SECS 

DATA  SAMPLING  PERIOD  -  . 10667E-O1  SECS 

POSITION  OF  THE  SIDESLIP  AND  INCIDENCE  ANGLE  SENSORS 

XBETA  «  . OOOOOE+OO  ZBETA  =  .OOOOOE+OO 

X ALPHA  -  . 13862E+02  YALPHA  -  .OOOOOE+OO 
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**  INITIAL  STATE  ESTIMATE  ** 

************************************************ 

*  VARIABLE  *  ESTIMATE  *  VARIANCE  * 

************************************************ 

*  U  *  . 12304E+03  *  . 40000E+01  * 

************************************************ 

*  W  ♦  . 14342E+02  *  .40000E+01  * 

************************************************ 

*  THET  *  . 11188E+00  *  .40000E-03  * 

************************************************ 

*  H  *  . 39265E+03  *  .40000E+01  * 

************************************************ 

**  INITIAL  ESTIMATE  OF  PARAMETERS  ♦♦ 

****************************************************************** 

*  VARIABLE  *  ESTIMATE  *  VARIANCE  *  NOISE  VARIANCE  * 

****************************************************************** 

*  BAI  *  . OOOOOE+OO  *  . 60000E-01  *  .OOOOOE+OO  * 

****************************************************************** 

*  BAZ  *  .OOOOOE+OO  *  .BOOOOE-Ol  *  .OOOOOE+OO  * 

****************************************************************** 

♦BP  *  .OOOOOE+OO  ♦  . 50000E-03  ♦  .OOOOOE+OO  * 

****************************************************************** 

*  BVEL  *  .OOOOOE+OO  ♦  . lOOOOE+Ol  ♦  .OOOOOE+OO  * 

****************************************************************** 

*  BALP  ♦  .OOOOOE+OO  ♦  .60000E-03  ♦  .OOOOOE+OO  ♦ 

****************************************************************** 

*  BTHE  ♦  .OOOOOE+OO  ♦  .B0000E-03  ♦  .OOOOOE+OO  * 

****************************************************************** 

»  BH  *  .OOOOOE+OO  ♦  . lOOOOE+Ol  ♦  .OOOOOE+OO  * 

****************************************************************** 

*  LALP  ♦  .OOOOOE+OO  ♦  .50000E-03  ♦  .OOOOOE+OO  ♦ 

****************************************************************** 
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**  DATA  PROCESSING  NO.  5  ** 


**  PARAMETERS  ESTIMATE  ** 

******* ***************************************** 

*  *  FINAL  PREDICTION  BY  THE  EKF  * 

*  VARIABLE  ************************************* 

*  *  ESTIMATE  *  STANDARD  ERROR  * 

************************************************ 

*  BAX  *  - . 14460E+01  *  .21644E-01  * 

************************************************ 

*  BAZ  *  - . 875B7E-01  *  .41061E-02  * 

************************************************ 

*  BP  *  - . 84207E-03  *  .29816E-04  * 

************************************************ 

*  BVEL  *  - . 13971E+02  *  .21119E+00  * 

************************************************ 

*  BALP  *  . 93545E-01  *  . 87074E-03  * 

************************************************ 

*  BTHE  *  .74014E-01  *  .947B9E-03  * 

************************************************ 

*  BH  *  . 27258E+01  *  .41214E+00  » 

************************************************ 

*  LALP  *  -.20941E-01  *  99603E-03  * 

************************************************ 


Appendix  G  -  A  priori  Data  from  Six  Degree  of  Free 
dom  Flight  Dynamic  Model 


The  six  degree  of  freedom  flight  dynamic  mathematical  model  of  the  F-lllC  is  used 
to  obtain  a  priori  values  of  the  aerodynamic  stability  and  control  derivatives  for 
use  in  the  parameter  estimation  techniques.  This  program  stores  the  aerodynamic 
data  in  data  tables  in  both  derivative  and  coefficent  form.  The  coefficent  form 
is  employed  where  the  aerodynamic  forces  are  non-linear.  Additional  subroutines 
were  developed  to  extract  local  derivatives  from  the  non-linear  model  using  the 
ACSL  simulation  language  Jacobian  analysis  facility.  Details  of  this  procedure  are 
given  in  References  [21]  and  [22].  Configuration  data  and  initial  conditions  for  the 
simulation  are  contained  in  file  TABLEE.NEW  shown  below.  The  aircraft's  mass 
characteristics  obtained  from  program  CGCALC-P10R3  described  in  Appendix  B, 
are  presented  in  the  first  part  of  the  table. 


WEIGHTS 

!  AND  INERTIAS 

6 

QUANTITIES 

#(I) 

QUANTITY 

VALUE(F) 

#(I) 

QUANTITY 

VALUE (F) 

1 

WEIGHT 

69900.0000000 

2 

IYY 

344866 . 0000000 

3 

IXX 

73343.0000000 

4 

IZZ 

411017.0000000 

6 

IXZ 

376.0000000 

6 

BLANK 

. 0000000 

INITIAL  CONDITIONS  (ALFW  IS  OVERWRITTEN  BY  TRIM  ROUTINE) 

18 

QUANTITIES 

#(I)  QUANTITY 

VALUE (F) 

#(I) 

QUANTITY 

VALUE (F) 

1 

GAMMAO 

.0000000 

2 

PO 

. 0000000 

3 

ALFW 

3 . 6007390 

4 

QO 

.0000000 

6 

BETAO 

.0000000 

6 

RO 

. 0000000 

7 

TRACKO 

. 0000000 

8 

XO 

.0000000 

9 

PHIDO 

. 0000000 

10 

YO 

. 0000000 

11 

XMACHO 

.3160000 

12 

ALTO 

1000.0000000 

13 

VWN 

. 0000000 

14 

VWE 

. 0000000 

16 

VWD 

. 0000000 

16 

DGAMMA 

.0000000 

17 

DVR 

. 0000000 

18 

CDWPS 

. 0000000 

AIRCRAFT  GEOMETRY 

4 

QUANTITIES 

#(I)  QUANTITY 

VALUE (F) 

#(I) 

QUANTITY 

VALUE (F) 

1 

SWEEP 

16.0000000 

2 

XONC 

. 2446800 

3 

ZONC 

.0380852 

4 

BLANK 

. 0000000 
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With  file  TABLEE.NE  W  altered  to  suit  the  particular  case,  running  of  the  model 
can  commence  (>  indicates  user  response  ): 

>  fllllatlong 

Information  describing  the  aerodynamic  data  base  used  in  the  simulation  is 
output  to  the  screen  and  the  ACSL  commands  and  flags  required  to  be  set  by  the 
user  are  provided.  These  included  the  commands  required  to  set  up  the  initial  run 
time  parameters  and,  for  the  example  shown,  to  set  the  flap  deflection  to  25  degrees, 
the  trim  flag  to  true  to  ensure  the  equations  are  trimmed  in  equilibrium,  and  the 
derivative  finder  option  to  true. 

>  S  CMD  -  10 

>  S  FLAPS  -  26. 

>  S  CTTRIM  -  .T. 

>  S  DEFIND  -  . T. 


>  START 

The  current  status  of  the  main  switches  and  flags  defining  the  control  system 
settings,  the  analysis  facilities,  and  the  configuration  is  presented  below:  The  CON¬ 
VERGENCE  signal  indicates  that  trimming  of  the  aircraft  has  taken  place. 


***** ****** ******** ********************************** 

♦  CURRENT  STATUS  OF  SET  UP  SWITCHES  * 

***************************************************** 

*  CONTROLS  SYSTEM  SWITCHES  * 

***************************************************** 


* 

RDAMP 

*  ROLL  DAMPER 

* 

F 

* 

PDAMP 

*  PITCH  DAMPER 

* 

F 

* 

YD  AMP 

*  YAW  DAMPER 

* 

F 

* 

RCAUG 

*  ROLL  CONTROL  AUGMENTATION 

* 

F 

* 

PC  AUG 

*  PITCH  CONTROL  AUGMENTATION 

* 

F 

* 

STA 

*  SERIES  TRIM  ACTUATOR 

* 

F 

* 

* 

* 

* 

* 

* 


***************************************************** 
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***************************************************** 

*  GENERAL  SET  UP  SWITCHES  * 

***************************************************** 

*  RADAPT  *  ROLL  ADAPTIVE  GAIN  FINDER  *  F  * 

*  PADAPT  *  PITCH  ADAPTIVE  GAIN  FINDER  *  F  * 

*  DEFIND  *  AERODYNAMIC  DERIVATIVE  FINDER  *  T  * 

*  FLTDAT  *  FLIGHT  DATA  INPUT  *  F  * 

*  CTTRIM  *  TRIM  AIRCRAFT  WITH  THRUST  *  T  * 

*  PTRIM  *  TRIM  AIRCRAFT  WITH  ELEVATORS  *  F  * 

*  RTRIM  *  TRIM  WITH  STEADY  SIDESLIP  *  F  * 

*  DUMP  *  GENERATE  DEBUG  OUTPUT  *  F  * 

*  PLOTFL  *  PRODUCE  3-D  GRAPICS  FILE  *  F  * 

***************************************************** 

*  IMPORTANT  CONSTANTS  * 

***************************************************** 


*  XCG  *  CENTRE  OF  GRAVITY  POSITION  *  .24468  * 

*  RFB  *  ROLL  CONTROL  FEEDBACK  *  1.00000  * 

*  AYFB  *  LATERAL  ACCELEROMETER  F/BACK  *  1.00000  * 

***************************************************** 

*  TAKE  OFF  AND  LANDING  CONFIGURATION  * 

***************************************************** 

*  LANDGR  *  LANDING  GEAR  DOWN  *  F  * 

*  SPEEDB  *  SPEED  BRAKE  DOWN  *  F  * 

*  SLATS  *  LEADING  EDGE  SLATS  EXTENDED  *  T  * 

*  GRDEFF  *  IN  GROUND  EFFECT  *  F  * 

*  FLAPS  *  TRAILING  EDGE  FLAPS  EXTENSION  *25.00000  * 

***************************************************** 

CONVERGENCE 

>  STOP 


File  F111DERV.0UT  shown  below  contains  the  a  priori  derivatives  calculated 
by  the  simulation  program  in  a  format  suitable  for  use  in  the  input  file  FAR01.DAT 
used  by  program  MMLE3. 

*********NON-DIMENSIONAL  LONGITUDINAL  AERODYNAMIC  DERIVATIVES.  ********* 


CTv 

a 

-.2937838 

CLv 

a 

-.0001830 

CDv 

m 

-.2937219 

CMv 

m 

- .0106612 
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CLalpha  3 
CDalpha  = 
CMalpha  ■ 
CLdalpha  - 
CMdalpha  - 
CLq  ■ 

CMq 

*********non-dimensional 


.1316456 
. 3000140 
-.0101156 
2 . 5420000 
-4.6251057 
5.7770060 
-14.2470031 

LONGITUDINAL  CONTROL 


DERIVATIVES .  ********* 


CLeta  -  .0174604 
CData  -  - . 0020467 
CMeta  -  - . 0353430 


*»****»**NON-DIMENSIONAL  LATERAL  AERODYNAMIC  DERIVATIVES.  ********* 


CYbeta 

M 

-1.0761711 

Clbeta 

3 

-.0471412 

CNbeta 

S 

.0671741 

CYp 

m 

-.0606060 

Clp 

■ 

-.4218507 

CNp 

■ 

.0151035 

CYr 

m 

.3232533 

Clr 

m 

■0426302 

CNr 

m 

-.1817570 

CYdbeta 

m 

-.1101400 

Cldbeta 

3 

-.0100330 

CNdbeta 

= 

.0351486 

*********NON-DIMENSIONAL  LATERAL  CONTROL  DERIVATIVES.  ********* 


CYalleroa" 

.0726516 

Claileron- 

- . 0346603 

CNaileroa- 

- . 0307724 

CYrudder  ■ 

. 0084668 

Clrudder  ■ 

. 0007033 

CNrudder  * 

-.0227010 

CYspoiler* 

. 0046761 

Clapoiler* 

-.0386813 

CNspoiler* 

-.0004143 
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Appendix  H  -  Parameter  Identification  Program  MMLE3 

The  data  processing  described  in  previous  Appendices  is  carried  out  to  prepare  the 
flight  data  for  use  by  the  parameter  identification  program  MMLE3.  Considerable 
background  knowledge  is  needed  to  apply  the  program.  Only  the  program  operation 
procedures  will  be  covered  in  this  Appendix.  A  detailed  description  of  the  program 
MMLE3  is  given  in  Reference  (15].  Shown  below  is  an  example  of  an  input  file 
PAR01.DAT  which  provides  all  the  data  required  for  the  third  run  of  a  phase  2 
flight,  event  72,  along  with  a  brief  comment. 

STANDARD  AIRCRAFT  ROUTINES 

END  ONCE 

NEW 

H1620036.0V3  18  TS1601  MMLE3  F111C  LONG  STAB  P2F1E72  RUN  3  SLATS  36 

FLAPS 
$WIND, 

AREA-660 . CH0RD=8 . 8 , SPAN-70 . 0 , 

CG-. 24609, 

STAB-T, 

PRINT  =F . 


Angle  of  attack  scale  factor  Ka  and  instrument  locations 


KALF-1 .06, 


XALF 

a 

46.6107, 

YALF 

=  .0000 

ZALF 

= 

2.6636 

XB 

3 

44.7032, 

YB 

=>  . 0000 

ZB 

= 

3.0618 

XAX 

- 

2.6241, 

YAX 

-  -2.4308 

ZAX 

- 

1.3877 

XAY 

3 

2.4218, 

YAY 

-  -2.4683 

ZAY 

= 

1.3877 

XAN 

3 

20.8674, 

YAN 

-  .0000 

ZAN 

3 

1.6610 

SEND 

0. 

0. 

0. 


A  priori  aerodynamic  information  including  matrix  element  positions  and  values 
The  default  longitudinal  AN  matrix  is 

'  cLa  c£f  0  ■ 

Cma  Cm,  0 
0  0  0 

The  default  lateral  AN  matrix  is 
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CY? 

°y t 

CYr 

0 

Ci, 

Q, 

0 

Cup 

Cn, 

Cn, 

0 

0 

0 

0 

0 

The  default  longitudinal  BN  matrix  is 


f  c<-.. 


The  default  lateral  BN  matrix  is 

'  CV4r 

ctf,  clh 

Cnia  C„Sr 
0  0 


CMA  LONG  AN(2, 1) 

-.01012 

CMQ  LONG  AN(2,2) 

-14.24700 

CMAD  LONG  RN(2,1) 

-4.62611 

CMDE  LONG  BN (2.1) 

- . 03634 

CNA  LONG  CN(4 . 1) 

. 13166 

CNQ  LONG  CN(4.2) 

6.77710 

CNDE  LONG  DN(4,1) 

.01746 

END 

H1620036.0V3  16  TS1601  MMLE3  F111C  LONG  STAB  P3F1E72  RUN  3  SLATS  36 

FLAPS 

$USER  LONG  -T. 

KALF-1.06 

W-  60000.86030  ,CG-  .24468, 

IX-  73343.67180  ,IY-  344866.78100, 

IZ-41 1017. 06200  ,IXZ-  3076.78710, 


XALF  - 

46.6107, 

YALF  - 

.0000  , 

ZALF  - 

2.6635 

XB 

44.7032, 

YB  - 

.0000  , 

ZB  - 

3.0618 

XAX  - 

2.6241, 

YAX  - 

-2.4308  , 

ZAX  - 

1 . 3877 
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XAY  =  2.4216,  YAY  =  -2.4683  .  ZAY 

XAN  -  20.6674,  YAN  =■  .0000  ,  ZAN 

SEND 


1.3877 

1.6610 


Namelist  INPUT  is  described  in  detail  in  Reference  [15]  section  3.3.8  with  most 
of  the  variables  being  equal  to  a  real  or  integer  type  number. 

$ INPUT  CARD-T , NREC=24 , ERRTH-F . 

NUPLT-1 .NEXPLT-6 , 

ERRMAX-10.E+30 

TIMVAR-T, 

SPS-60.TEST-F. 

MZ-4, 

MU-1, 

MX- 3 , 

NQITER-29 , 

ITG-30, 

WAPR-1.0, 

ITAPR-20, 

ZSCALE  (6) -67 . 296 
ZSCALE(1)-1 .0 
SEND 

100000000  100006000 

Some  of  the  aerodynamic  parameters  are  linked  through  hard  constraints  to  re¬ 
duce  the  number  of  independent  variables  and  to  improve  the  identification  process. 
For  the  longitudinal  motion  the  following  constraints  based  on  the  a  priori  model 
values  were  used: 

Cmi  —  0.32464  x  cmj 
Q,  =  330.86868  x  C„4< 

Iio,  ~ 

CLa  ~  Csa 


HARD  5 

HARD  6 

RN(2 , 1)  -  AN(2,2) 

*  . 32464 

AN (1,2)  -  DN(4 , 1) 

*  330.86868 

CN(4,2)  -  DN(4 , 1) 

*  330.86868 

BN(1 , 1)  -  DN(4 , 1) 

*1.0 
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AN(l.l)  =  CM (4,1) 
GGI  6 

10.  80. 
END 


Flight  data  is  provided  in  the  following  order: 
time,  a,q,V,9,an,q,ax,Se,Sc,  Silon>,tud,nal >  S2,aat„udinai >  ^ q,0,p,r,ay,p,r, 

Sat  Sr  l  8  llal.ral  >  ^2|ai.raJ 


100000000 

-.0336 

2.6800 

.0948 

366.6670 

1.1361 

1 . 0669 

1 . 6889 

-4.4976 

146.0993 

.0000 

.0000 

.0000 

-1.6664 

1198.0884 

.3209 

.2061 

-.5241 

-.6104 

-1.1277 

-.6927 

-.0210 

1 . 4067 

-10.7161 

.8083 

100000017 

-.0331 

2.0632 

.0423 

356.2252 

1.1361 

1.0628 

9.8166 

-4.6087 

145.7467 

.0000 

.0000 

.0000 

-1.6664 

1197.8882 

.3206 

-.0678 

-.4116 

-.6326 

-1.1277 

- .  6927 

- . 0247 

9.6135 

-6.4647 

.7972 

100000034 

- .0346 

1.9216 

.0948 

356.6670 

1.2242 

1.0161 

3.7826 

The  results  of  the  identification  process  are  written  into  file  FAR03.DAT  and 
shown  below  are  the  results  obtained  after  12  iterations: 


*1.0 


200. 


0.0 


100. 


WEIGHTED  ERRORS 
ERROR  SUM  *  3.0000 

ALPHA  q  THETA  AN 

1.000  1.000  .0000  1.000 


DET  (RSq)  -  -3.404288 

AN  3 

.1202  6.500 

- . 1918E-01  -27.80 

. 0000  . 0000 


BY  3 
.0000 
.0000 
.0000 
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t 


i 

f 

I 

t 

I 

► 


BN 

. 2190E-01 

3 

BY 

1 

- . 3046E-01 

.0000 

SN 

.7031 

3 

BY 

1 

- .  7449E-01 

.0000 

.0000 

CN 

4 

BY 

3 

1.160 

62.34 

.0000 

.0000 

1.000 

.0000 

.0000 

.0000 

1.000 

.r?02 

6.609 

.0000 

DN 

.0000 

.0000 

.0000 

.2190E-01 

4 

BY 

1 

HN 

4 

BY 

1 

.0000 

.0000 

.0000 

.7113 

WEIGHTED  ERRORS 
ERROR  SUM  -  2.9084 

ALPHA  q  THETA  AN 

1.008  1.021  .0000  .9691 

DET  (RSq)  »  -3.485601 

GGI  4  BY  4 


4.685 

.0000 

.0000 

.0000 

.0000 

1.185 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

434.5 

ITERATION  12  COMPLETED.  H1620035.0V3  16  TS1691  MMLE3  F111C  LONG  STAB 

P3F1E72  RUN  3  SLATS  35  FLAPS 

COST  FUNCTION  CONVERGED  WITHIN  . 10E-02  BOUND. 


WEIGHTED  ERRORS 
ERROR  SUM  -  3.0000 
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ALPHA 

1.000 


Q 

1.000 


THETA 

.0000 


AN 

1.000 


\ 


•  '  V 


/ 


DET  (RSq)  -  -3.485601 


CRAMER -RAO  BOUNDS. 

AC 

3  BY 

3 

.0000 

.0000 

.0000 

. 1363E-03 

.7636 

.0000 

.0000 

.0000 

.0000 

BC 

.0000 

. 3763E-03 
.0000 

3  BY 

1 

SC 

. 8623E-02 
. 197 IE-02 
.0000 

3  BY 

1 

CC 

4  BY 

3 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

DC 

. 8800E-03 

.0000 

,0000 

.0000 

.0000 

.0000 
.  8664E-03 

4  BY 

1 

HC 

.0000 

.0000 

.0000 

■ 6932E-02 

4  BY 
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u  ■ » 


T*- 

s 


w 


Appendix  I  -  Longitudinal  Equations  of  Motion 

State  variables  x  =  (a ,q,9) 

Control  variables  u  =  (Se) 

Observation  variables  r  =  (am>qm,9m,anm,aXm]qm) 

The  nonlinear  longitudinal  state  equations  are  : 

d  =  +  do)  +  q  +  ^-R(coa9coaa  +  sinasinfl) 

mV  V 

T  n  . 

- —  R  sin  ct 

mV 

lvyq  =  qScRCm 
9  =  q  +  0o 

The  do  and  9q  are  included  to  allow  for  instrument  biases. 


The  longitudinal  observation  equations  are  : 

<*m  =  Ka{a-~-q) 


qm  —  q 


=  9 


qS  r*,  x<3n  •  ,  za.. 

mg  N  gR  R2g 


.,2 


qS 


"  ~^rgCA  +  jRq~Wgq  +TTg 

9m  =  9  +  90 


The  qo  is  the  instrument  bias  on  q  . 
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The  expansions  of  the  longitudinal  force  and  moment  coefficients  are  : 
Cn  =  CjVaa  -r  -r  Cs6Jt  -r  C/v0 


<70  ac 

Cm  =  Cm<1a  r  Cm?  ^  f  Cmj>  6e  +  Cmo  +  Cm<i  ^7^ 


CA  -  CAaa -r  Cyi,7rrr -i- Ca^^  +  Cx0 


Cl  =  C^y  cos  a  -  Ca  sin  a 
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7 


r 

j 

i 

Appendix  J  -  Lateral  Equations  of  Motion 

•  State  variables  x  —  {P,p,r,p) 

a 

Control  variables  u  =  (<5a,  <^r>  &tp) 

Observation  variables  2  =  (/ 3m,pm,rm,4>m,aVm ) 

The  nonlinear  lateral-directional  state  equations  are  : 

qS  g 

=  — —R(Cy  +  A))  +  77  cos  0  sin  0  +  p  sin  a  -  r  cos  a 
mV  V 

plz  -  rlzz  =  qSbRCt^^(Iv-  I,). rW| 

r/.-p/«,  =  qSbRCn  -r  ^{Iz  -  /„)  - 

<p  =  p  +  r  cos  <p  tan  0  +  g  sin  <p  tan  0  + 

The  terms  /?q  and  0o  are  included  to  allow  for  instrument  biases. 


The  lateral  observation  equations  are  : 
Pm 


Pm 

=  p 

rm 

=  r 

0m 

=  4> 

aVm 

=  ficy- 

mg 

Pm 

=  P  +  P0 

^  m 

—  r  +  tq 

gRV  gR  R*g 


(P2  +  r2) 


The  terms  po  and  fo  are  the  instrument  biases  on  p  and  r. 


The  expansions  of  the  lateral  force  and  moment  coefficients  are  : 

cy  = 

Cl  = 

Cn  = 

where  the  S  term  is  summed  over  all  controls.  For  the  50°  and  72.5°  sweeps,  the 
aileron  (6a)  and  rudder  (6r)  are  used.  For  sweeps  less  than  or  equal  to  47°,  the 
spoilers  (6,p)  are  also  operative. 


Cy,P  +  CY, ^  +  Cy-  2YR  +  Cy‘S  +  Ck° 


pb  rb  3b 

CipP  +  C‘rzy~R  +  Ct'2VR  +  C‘sS  +  C'°  +  C|*; 


t  2  VR 


pb  rb  3b 

Cnff3  +  CnP2VR  +  Cn'2VR  +Cn*S  +  Cno  +  CnHw~R 


2 VR 
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Appendix  K  -  Output  File  Notation 


The  notation  for  naming  and  cataloging  MMLE3  results  uses  the  following  infor¬ 


mation:  case  type,  sweep,  altitude,  Mach  numbet  and  version. 


A 

casetype 


sweep 

50 


£30 

altitude 


Machnumber 

M  1.2 


V2 


version 


Case  type 


A  Longitudinal  MMLE3 
B  Lateral  MMLE3 
C  Longitudinal  MMLE3P 
D  Lateral  MMLE3P 
E  Longitudinal  MODEL 
F  Lateral  MODEL 


Where  MMLE3  is  the  standard  parameter  estimation  program.  MMLE3P  is  a 
version  with  fixed  weightings  and  MODEL  contains  the  a  priori  results  obtained 
from  the  six  degree  of  freedom  flight  dynamic  model  of  the  F-111C. 


sweep 

16 

16  degrees 

26 

26  degrees 

35 

35  degrees 

45 

45  degrees 

50 

50  degrees 

72 

72.5  degree 

Height 

H05 

5000  feet 

H10 

10000  feet 

H20 

20000  feet 

H30 

30000  feet 

H40 

40000  feet 

H50 

50000  feet 
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Version  longitudinal 


1 

Phase  1  +ve  g  pullup 

2 

Phase  1  -ve  g  pushover 

3 

Phase  2  +ve  g  pullup 

4 

Phase  2  -ve  g  pushover 

5 

Phase  1  -t-ve  g  pullup  repeat 

6 

Phase  1  -ve  g  pushover  repeat 

7 

Phase  2  +ve  g  pullup  repeat 

8 

Phase  2  -ve  g  pushover  repeat 

lateral 

1 

Phase  1  initial  roll  to  the  right 

2 

Phase  1  initial  roll  to  the  left 

3 

Phase  2  initial  roll  to  the  right 

4 

Phase  2  initial  roll  to  the  left 

5 

Phase  1  initial  roll  to  the  right  repeat 

6 

Phase  1  initial  roil  to  the  left  repeat 

7 

Phase  2  initial  roll  to  the  right  repeat 

8 

Phase  2  initial  roll  to  the  left  repeat 

Longitudinal  Analysis 

Shown  below  is  a  typical  longitudinal  output  file  produced  by  MMLE3  with  the 
following  information: 

File  designation,  title,  flag  to  indicate  Cm„  has  been  corrected  for  C  of  G  position, 
actual  sweep  angle,  actual  Mach  number,  actual  altitude,  actual  all  up  weight  in 
pounds,  C  of  G  as  a  %  of  mean  aerodynamic  chord. 


A60H30M1 . 2V2  50  TS1691  MMLE3  F111C  LONG  STAB  L8E01  RUN  8_NG 

CMACORR  50  +1.2287  +30076.6113  +74812.6250  +0.2978 

Estimates  for  the  longitudinal  stability  and  control  derivatives  are  presented 
followed  by  their  Cramer-Rao  bounds. 


Cl, 

Cma 

Cm, 

Cm<j 

Cls. 

Cm« 

+0.0747  +6.7827 

-0.0565  -17.0710 


68 


-3.2007 

+0.0107  -0.0264 

+0 . 0007  +0 . 0000 

+0.0001  +0.2888 

+0 . 0000 

+0.0004  +0.0001 


The  last  line  contains  the  weightings  for  the  output  variables  a,q,9,an  used  in 
the  cost  function  equation. 


+23.2331  +16.0489  +0.0000  +1081.3801 


Lateral  Analysis 

Shown  below  is  a  typical  lateral  output  file  produced  by  MMLE3  with  the  following 
information: 

File  designation,  title,  actual  sweep  angle,  actual  Mach  number,  actual  altitude, 
actual  all  up  weight  in  pounds. 


D72h20ml.6v2  72.6  plflle4a 

72.5  +1.62387714  +20636.6234  +72698.6 


Estimates  for  the  lateral  stability  and  control  derivatives  are  presented  followed 
by  their  Cramer-Rao  bounds. 


Cy, 

Cy, 

Cy, 

Cb 

Ci, 

Crip 

Cn, 

Cn, 

cy <. 

C*t  r 

CY* 

V'r 

Cl* 

Cnlt 

Cnfr 

Cn*  . 

-8 . 48778623E-03 
-6.61818797E-04 
+0 . 76631927E-04 
+1 . 68686606E-03 
-8.60938463E-04 
-1 . 13606608E-04 
+1.06449871E-04 
+8 . 6396003E-06 


-3 . 59000004E-02 
-4.74088229E-02 
+2 . 00000009E-03 
+7 . 39969196E-04 
+4.41243264E-06 
-4 . 27318164E-04 
+0.0 

+4 . 21166204E-04 


+0.286699993 

+2.06999996E-O2 

-0.184223116 

+0.0 

+0.0 

+0.0 

+0.0 

+0.0 
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+5 . 87599834E-06 
+4 . 32894968E-06 
+3.71340024E-06 
+2.76041328E-06 


+0.0 

+8.66981717E-05 
+6.69186238E-06 
+5 . 60222906E-08 


+2.9275713E-03 

+0.0 

+0.0 

+0.0 


The  last  line  contains  the  weightings  for  the  output  variables  l3,p,r,<p  and  av  in 
the  cost  function  equation. 


+6.00899982  +0.400000005  +15.0 


+0.0 


+16000.0 


Take-Off  and  Landing  Analysis 


The  notation  for  the  take-off  and  landing  configuration  is  as  follows: 

sweep  flapsetting 


H  16 

casetype 


240 

VI  AS  knots 


VI 

version 


Case  type  H  Longitudinal  MMLE3 

I  Lateral  MMLE3 

J  Longitudinal  MMLE3P 

K  Lateral  MMLE3P 

L  Longitudinal  MODEL 

M  Lateral  MODEL 

VI AS  150  150  knots  indicated 

160  160  Knots  indicated 

180  180  Knots  indicated 

200  200  Knots  indicated 

215  215  Knots  indicated 

240  240  Knots  indicated 

Flap  setting  00  0  flaps,  slats  extended 

15  15  degrees  flap,  slats  extended 

25  25  degrees  flap,  slats  extended 

35  35  degrees  flap,  slats  extended 

Version  designation  is  identical  for  the  clean  aircraft  configuration  and  the  take 
off  and  landing  configuration. 


Table  1:  Instrumentation  Channels  used  for  Flight  Dynamic  Analysis 


Quantity 

Symbol 

Units 

Range 

Accuracy 

ACCEL  LONG  CG* 

a* 

g 

±5 

±0.05 

ACCEL  LAT  CG 

av 

g 

±5 

±0.05 

ACCEL  VERT  CG** 

a* 

g 

±10 

±0.05 

ROLL  RATE 

P 

deg/ 8  ec 

±300 

±3 

PITCH  RATE 

9 

deg/sec 

±100 

±i 

YAW  RATE 

r 

deg/  a  ec 

±50 

±0.5 

ROLL  ACC  CG 

P 

rad/  sec2 

±10 

±0.05 

PITCH  ACC  CG 

7 

rad/sec2 

±5 

±0.05 

YAW  ACC  CG 

r 

rad/aec 2 

±5 

±0.05 

ROLL  ANGLE 

4> 

deg. 

±180 

±0.5 

PITCH  ANGLE 

e 

deg. 

±180 

±0.5 

ANGLE  OF  ATTACK* 

a 

deg. 

-3  —  25 

±0.5 

ANGLE  OF  SIDESLIP* 

0 

deg. 

±24 

±0.5 

VELOCITY 

V 

kt. 

0  —  900 

±10 

MACH  No.* 

M 

- 

0.3  —  1.8 

±0.001 

ALTITUDE* 

H 

ft. 

-500  —  55000 

±1.5 

WING  SWEEP 

A 

deg. 

16  -  72.5 

±0.05 

STABILATOR  (right) 

deg. 

-30  —  15 

±0.1 

STABILATOR  (left) 

deg. 

-30  —  15 

±0.1 

RUDDER 

Sr 

deg. 

±30 

±0.1 

SPOILER  (right) 

deg. 

0  —  45 

±0.1 

SPOILER  (left) 

deg. 

0  —  45 

±0.1 

STICK  POS  (long) 

in. 

-4.4  —  3.6 

±0.05 

STICK  POS  (lat) 

in. 

±5 

±0.05 

RUDDER  PED  POS 

in. 

±3 

±0.03 

*  NBTU  measurements  available  in  Phase  2 
t  Altitude  and  Mach  no.  coarse  and  fine  readings  available 
|  Crew  module  normal  accelerometer  used  in  Sights  2  and  3 
of  Phase  2  because  of  a  fault  in  the  CG  accelerometer  signal 


71 


Table  3:  Take-Off  and  Landing  Aircraft  Configurations 


Table  4:  Supplementary  Manoeuvres 


Manoeuvre 


Phase 


Altitude 


Mach  no./ 


Sweep 


Roller  Coaster 


1 


Lateral  Oscillations 


2 

1 


Dutch  Roll 


2 


Steady  Heading  Sideslip 
Right 

Left _ 

Steady  Level  Trim 


(ft.) 

4000 

5000 

10000 

10000 

1000 

2000 

5000 

10000 

20000 

10000 


20000 


KCAS 
0.60 
0.60 
0.80 
0.80 
1.20 
0.70 
0.40 
143  kn. 
0.60 
0.80 
0.80 
1.20 
0.60 
0.70 
0.70 
0.70 
0.65 
0.74 


(deg.) 

30 

30 

56 

55 
72.5 

35 

16 

16 

30 

56 
55 
54 
26 
26 
35 
45 
26 
26 


0.60 

0.60 

0.60 


2 

2 

2 


10000 

10000 

10000 


26 

26 

26 


Table  5:  Weighing  Information  for  Phase  1  Aircraft  ,  16°  Sweep  Flaps 
Down 


Configuration 


Wing  Sweep:  16  Degrees 
Flaps:  Fully  extended 

Slats:  Fully  extended 

Gear  :  Down 


Support 

Points 

Weight 

(lb) 

Port 

Starboard 

Nose 

22208 

22357 

6108 

Support  Point 
Total 

Weight 

(lb) 

Arm 

(in) 

Moment 
(lb- in) 

Main 

44565 

561.4 

25018791 

Nose 

6108 

264.8 

1617398 

Total 

50673 

525.6 

26636189 

Table  6:  Weighing  Information  for  Phase  1  Aircraft  ,  16°  Sweep 


Configuration 

Wing  Sweep:  16  Degrees 
Flaps:  Up 

Slats:  Up 

Gear  :  Down 


Support 

Points 

Weight 

(lb) 

Port 

Starboard 

Nose 

22207 

22242 

6158 

Support  Point 
Total 

Weight 

(lb) 

Arm 

(in) 

Moment 

(lb-in) 

Main 

44449 

561.4 

24953669 

Nose 

6158 

264.8 

1630638 

Total 

50607 

525.3 

26584307 
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Table  7:  Weighing  Information  for  Phase  1  Aircraft  ,  26°  Sweep 


Configuration 

Wing  Sweep:  26  Degrees 
Flaps:  Up 

Slats:  Up 

Gear  :  Down 


Support  Point 
Total 

Weight 

(lb) 

Arm 

(in) 

Moment 

(lb-in) 

MaiD 

44984 

561.4 

25254018 

Nose 

5713 

264.8 

1512802 

Total 

50697 

528.0 

26766820 

Support 

Points 

Weight 

(lb) 

Port 

Starboard 

Nose 

22445 

22539 

5713 

Table  8:  Weighing  Information  for  Phase  1  Aircraft  ,  35°  Sweep 


Configuration 

Wing  Sweep:  35  Degrees 
Flaps:  Up 

Slats:  Up 

Gear  :  Down 


Support  Point 
Total 

Weight 

(lb) 

Arm 

(in) 

Moment 

(lb-in) 

Main 

45407 

561.4 

25491489 

Nose 

5233 

264  3 

1385698 

Total 

50640 

530.8 

26877188 

Support 

Points 

Weight 

(lb) 

Port 

Starboard 

Nose 

22631 

22776 

5233 
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Table  9:  Weighing  Information  for  Phase  1  Aircraft  ,  50°  Sweep 


Configuration 

Wing  Sweep:  50  Degrees 
Flaps:  Up 

Slats:  Up 

Gear  :  Down 


Support 

Points 

Weight 

(lb) 

Port 

Starboard 

Nose 

22631 

22776 

5233 

Support  Point 
Total 

Weight 

(lb) 

Arm 

(in) 

Moment 
(lb- in) 

Main 

45407 

561.4 

25491489 

Nose 

5233 

264.8 

1385698 

Total 

50640 

530.8 

26877183 

Table  10:  Weighing  Information  for  Phase  2  Aircraft  ,  26°  Sweep 


Configuration 

Wing  Sweep:  26  Deg 
Flaps:  Up 

Slats:  Up 

Gear  :  Down 

NBTU  installed 


Support  Point 
Total 

Weight 

(lb) 

Arm 

(in) 

Moment 
(lb- in) 

Main 

44520 

562.64 

25048732 

Nose 

6555 

269.55 

1766900 

Misc. 

23 

16998 

Total 

51098 

524.8 

26817331 

Support 

Weight 

Points 

(lb) 

Port 

22199 

Starboard 

22321 

Nose 

6555 

Not?  : 

Misc.  refers  to  an  adjustment  that  has  to  be  made  to  the  final  weight  and  moment 
arm  of  the  aircraft  to  compensate  for 

•  Any  items  weighed  but  are  not  part  of  the  basic  configuration  . 

•  Items  that  are  part  of  the  basic  configuration  that  were  not  in  the  aircraft 
when  it  was  weighed. 
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Figure  1:  Instrumented  F-111C  Aircraft  (A8-132)  Operated  by  ARDU  for 
Flight  Trial 


Digital  processing 

ELECTRONICS  BOX 


j  Flight  data  on  AFTRAS  tape 

V _ ) 


Section  4.1 


Sections  4.2 
and  4.3 


Section  4.4 


Section  4.5 


Section  4.6 


Section  4..' 


Figure  3:  Summary  of  Flight  Data  Processing  and  Analysis  Procedures 


ELXSI  __  VAX  750 


INPUT  FILES 


PROGRAMS 


OUTPUT  FILES 


INDATA.DAT 

OUTPUT.DAT 


EXTRACT 


LISTPARM 


OUTPUT.DAT 

LISTPARM.DAT 


LISTPARM.DAT 

PLOTS.DAT 


r  1 

[  TRPLOT 
— 
TRANS 

V  > 


PLOTS.DAT 

PLOTS. PLO 


DATABASE  ;  cGCALC  P10R3 
DGTABLE  r  l'-/' 


LISTPARM.DAT  I  FDP 


DATIN.ALL 

DATIN.DAT 


FPRSORT 


DATIN.OUT 


FPR 

v _ y 


MMLE3.PARA 

F111C.DAT 
FPR. LON 
FPR.LAT 
MMLE3.DAT 

DATIN.OUT 

FLINP.DAT 

FPRTOMMLE3 

PREOUT 

FLI1A.DAT 

STATE.DAT 

COMPAT.DAT 


Figure  4:  Summary  of  Computer  Programs,  Input  and  Output  File  Names 
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ELXSI 


INPUT  FILES 


PROGRAMS 


OUTPUT  FILES 


STATE.DAT 

FLI1A.DAT 

COMPAT.DAT 


r 


FPRPLOT 


METADAT 


TABLEE.NEW 

SETUP 

F111LONGDB 

F111LATDB 

F111TOLADB 


F1 1 1 LATLONG 


F111DERV 
METADAT 
F1 11  LATLONG. L 


FAR01.DAT 


MMLE3 

v _ _  J 


FOR03.DAT 

METADAT 

'RESULTFILE’ 


CM 

CO 

CL 


’RESULTFILE’ 


Fill 


'RESULTFILE'. DAT 
’RESULTFILE’. GRF 
'RESULTFILE'. PRN 


Figure  5.  Summary  of  Computer  Programs,  Input  and  Output  File  Names 
(continued) 


83 


Control 

input 


r 


v 


Noise 


Test 

aircraft 


Measured 

response 

- >1 
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Figure  7:  Block  Diagram  of  Maximum  Likelihood  Estimation, Procedure 

\ 


85 


AERODYNAMIC  DERIVATIVES 

used  in  Linear  Identification  Models 


Longitudinal 


Lateral 


Static  Derivatives 


s. 


Dynamic  Derivatives 


s. 


Sr 


Control  Derivatives 


S, 


CY  Sf  S' 

5  5  9 


m 

s 


Totals 


19 


Importance  of  derivatives  to  aircraft  motion. 
I"!  Prime 

Secondary 

Least 


Figure  8:  Summary  of  Aerodynamic  Parameters  Used  in  Linear  Identifi¬ 
cation  Models 
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